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An analytical and nonperturbative approach to SU(2) and SU(3) Yang-Mills thermo- 
dynamics is developed and applied. Each theory comes in three phases: A deconfining, 
a preconfining, and a confining one. We show how macroscopic and inert scalar fields 
emerge in each phase and how they determine the ground-state physics and the proper- 
ties of the excitations. While the excitations in the deconfining and preconfining phase 
are massless or massive gauge modes of spin 1 the excitations in the confining phase 
are massless or massive spin-1/2 fermions. The nature of the two phase transitions is 
investigated for each theory. We compute the temperature evolution of thermodynami- 
cal quantities in the deconfining and preconfining phase, and we estimate the density of 
states in the confining phase. Some implications for particle physics and cosmology are 
discussed. 
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1. Introduction 

The beauty, richness and usefulness of nonabelian gauge theories is generally appre- 
ciated. Yet, in a perturbative approach to gauge theories like the Standard Model 
of particle physics (SM) and its (non)supersymmetric extensions it is hard if not 
impossible to convincingly address a number of recent and not so recent experimen- 
tal and observational results in particle physics and cosmology: Nondetection of the 
Higgs particle at LEP ™ indications for a rapid thermalization and strongly collec- 
tive behavior of the plasma that emerges in the early stage of an ultra-relativistic 

IOIQI 

heavy-ion collision despite the fact that the ideal hydrodynamical expansion es- 
sentially obeys a free-gas equation of state, dark energy and dark matter, a strongly 
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favored epoch of cosmological inflation | 4 | 7 | 5 | 6 | m ^ ne verv early Universe l ^ l ^ l ^ l 
and today's accelerated cosmological expansion, a puzzling large-angle signal in the 
power spectrum of the cosmic microwave background for the cross correlation of 
electric-field polarization and temperature fluctuation s J^j, t he likely existence of 
intergalactic magnetic fields of so far unclarified origin HSL L4J an c [ the departure of 
~ 3 x 10 43 protons p. a. (and hardly any negative charges) ^^3J from the sun's sur- 
face (solar wind) which clearly is contradicting the charge conservation inherent in 
the SM. An analytical and nonperturbative approach to strongly interacting gauge 
theories may further our understanding of these phenomena. 

The objective of the present work is the thermodynamics of SU(2) and SU(3) 
Yang-Mills theories in four dimensions. It is difficult to gain insights in the dynamics 
of a strongly interacting four-dimensional gauge theory by analytical means if this 
dynamics is not severely constrained by certain global symmetries. We conjecture 
with Ref. 1 ^ that a thcrmodynamical approach is an appropriate starting point for 
such an endeavor. On the one hand, this conjecture is reasonable since a strongly 
interacting system, being in equilibrium, communicates distortions almost instan- 
taneously by rigid correlations, and thus a return to equilibrium takes place very 
rapidly. On the other hand, it turns out that the requirement of thermalization al- 
lows for an analytical and nonperturbative derivation of macroscopic ground states 
and the properties of their (quasiparticle) excitations in two of the three phases 
of each theory. A breakdown of equilibrium in a transition to the third phase is 
unproblematic since the dynamics then is uniquely determined by the remaining 
symmetry. 

Let us very briefly recall some aspects of the analytical approaches to ther- 
mal SU(iV) Yang-M ills th eory as they are discussed in the literature. Because of 
asymptotic freedom | 16 | 17 | one wou i,^ naively expect thermal perturbation theory 
to work well for temperatures T much larger than the Yang-Mills scale A since 
the gauge coupling constant g(T) logarithmically approaches zero for ? — > oo. It 
is known for a long time that this expectation is too optimistic since at any tem- 
perature perturbation theory is plagued by i nsta bilities arising from the infrared 
sector (weakly screened, soft magnetic modes ^j). As a consequence, the pressure 
P can be computed perturbatively only up to (and including) order g 5 . The effects 
of resummations of one-loop diagrams (hard thermal loops), which rely on a scale 
separation furnished by the small value of the coupling constant g, are s um marized 
in terms of a nonlocal effective theory for soft and semi-hard modes In the 
computation of radiative corrections, based on this effective theory, infrared effects 
due to soft magnetic modes still appear in an uncontrolled manner. This has led to 
the construction of an effective theory where soft modes are collectively described 
in terms of class ical fie lds whose dynamics is influenced by integrated semi-hard 
and hard modes I 20 | 21 | j n Q uan t um Chromo dynamics (QCD) a perturbative cal- 
culation of P was pushed up to order g 6 log q, and an additive 'nonperturbative' 
term at this order was fitted to lattice results ^1 Within the perturbative orders 
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a poor convergence of the expansion i s ob served for temperatures not much larger 
than the MS scale. While the work in^lis a computational masterpiece it could, 
by definition, not shed light on the missing, nonperturbative dynamics of the in- 
frared sector. Screened perturbation theory, which relies on a split of the tree- level 
Yang-Mills action by the introduction of variational parameters, is a very interest- 
ing idea. Unfortunately, this approach generates temperature dependent ultraviolet 
divergences in its presently used form, see ^51 for a recent review. 

The purpose of the present work is to report, in a detailed way, on the devel- 
opment of a nonperturbative and analytical approach to the thermodynamics of 
SU(2) and SU(3) Yang-Mills theory (see for intermediate stages). The reason 
why we consider only these two gauge groups is that for SU(iV) with N > 4 it is 
likely that the phase structure of the theory is not unique: In contrast to SU(2) 
and SU(3), which possess one confining (center), one preconfining (magnetic), and 
one deconfining (electric) phase, more than three phases may exist for an SU(A) 
Yang-Mills theory with N > 4. 

Our starting point is the derivation of a macroscopic ground state in the de- 
confining phase. This ground state originates from instantaneous, long-range corre- 
lations mediated by field configurations of topological charge ±1. Technically, this 
situation is described by a spatially homogeneous, quantum mechanically and sta- 
tistically inert scalar field 4>, which transforms under the adjoint representation of 
the gauge group, and a pure-gauge configuration of trivial topology solving the 
Yang-Mills equations subject to a source term provided by 4>. Both the field <j> and 
the pure-gauge configuration emerge after a spatial coarse-graining over quantum 
fluctuations down to a resolution corresponding to the length scale |</>| _1 . While <fi 
represents the spatial average over BPS saturated, topological defects, that is, 'large' 
quantum fluctuations the pure-gauge configuration is a manifestation of averaged- 
over plane- wave quantum fluctuations. 

Conceptually, our approach is similar to the macros copic Landau-Ginzburg- 
Abrikosov theory for superconductivity in metals | 25 | 26 | R eca [] that this theory 
describes the existence of a condensate of Cooper pairs in terms of a nonvanishing 
expectation for a complex scalar field tp (local order parameter), which is charged 
under the electromagnetic gauge group U(l), and in terms of a pure-gauge con- 
figuration. A nonzero value of ip is enforced by a phenomenologically introduced 
potential V. As a consequence, coarse-grained U(l) gauge-field modes <5a M (pho- 
tons), which are deprived of the microscopic gauge-field fluctuations contributing 
to the formation of Cooper pairs and their subsequent condensation, acquire a mass. 
Microscopically, the generation of a photon mass can be visualized as a large se- 
quence of elastic scattering processes off the electrons residing within individual 
Cooper pairs in the condensate. At a given photon momentum this slows down the 
effective velocity of propagation in comparison to a propagation without a Cooper- 
pairs condensate, see Fig. [I] In the superconducting phase the U(l) gauge symmetry 
is spontaneously broken by the Cooper-pair condensate, and physical phenomena 
associated with this breakdown can be analysed in dependence of the parameters 
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(a) 



(b) 

Fig. 1. Photon propagation without (a) and within (b) a Cooper-pair condensate. 



entering an effective action and in dependence of external conditions such as a 
magnetic field and/or temperature. 

When pursuing the idea of a dynamically generated, macroscopic ground state 
in each of the phases of an SU(2) or SU(3) Yang-Mills theory it turns out that a 
situation similar to superconductivity holds. Moreover, in a thermalized SU(2) or 
SU(3) Yang-Mills theory one is in the comfortable position of being able to derive 
the dynamics of macroscopic scalar fields from first principles. Thus the stabilizing 
potentials for each scalar field are uniquely determined (up to an undetermined mass 
parameter - the scale of the Yang-Mills theory). Each (gauge invariant) potential is 
specified by a unique microscopic definition of the scalar field's phase (in a suitably 
chosen gauge) and by the assumption that a dynamically generate d, cons tant mass 
scale exists. This assumption is supported by perturbation theory EilHI where the 
specification of the running of the gauge coupling requires such a boundary con- 
dition. The microscopic definition of the scalar field's phase is an average over an 
(nonlocal) operator saturated by noninteracting Bogomolnyi-Prasad-Sommerfield 
(BPSjjHS saturated configurations of topological charge ±1. If a particular phase 
supports propagating gauge modes then, in a second step, interactions between 
the topological defects and microscopic radiative corrections are taken into account 
macroscopically by pure-gauge solutions to the equations of motion for gauge fields 
residing in the trivial-topology sector. The source term for these equations of motion 
is provided by the (inert) scalar fields. 

More specifically, we have shown in SSI that at large temperatures upon spatial 

coarse-graining an adjoint scalar field <j> emerges due to the spatial correlations 

fool 

mediated by trivial-holonomy calorons . (By large temperature we mean large as 
compared to the dynamically generated scale.) We discuss in Sec. 12.3.41 why the 
critical temperature Tp for the onset of </>'s existence should be comparable to the 
cutoff-scale for a field theoretic description in fo ur d imensions. Trivial-holonomy 
calorons are BPS saturated (or selfdual) solutions I^QI to the Euclidean Yang-Mills 
equations at finite temperature. (The time coordinate r is compactified on a circle, 
< t < A.) The topological charge of these configurations is integer. (Whenever we 
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speak of a topological soliton in this section its antisoliton is also meant.) It will turn 
out that only calorons with topological charge one contribute to the moduli-space 
average defining 0's phase. The reason is that for higher-charge calorons the larger 
number of dimensionful moduli does not admit the definition of a dime nsio nless 
entity without the introduction of an explicit temperature dependence SSI. The 
latter, however, ought to be absent because of a temperature independent weight 
on the classical level. To understand macroscopic results generated by microscopic 
interactions between calorons and between their constituents an investigati on o f the 
properties of quantum corrected nontrivial-holonomy calorons is necessary S2. 

For a given SU(iV) gauge-field configuration nontrivial holonomy refers to the 
following property: Evaluate the Polyakov-loop on this configuration at spatial in- 
finity and observe that the result is not an clement 



1 exp 
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k 

N 



(fc = 0,---,JV-l) (1) 



of the center Zn of the gauge group. If this is the case then a mass scale exists in the 
configuration which determines the behavior Ai(|a:| — > oo). For a classical solution 
at finite temperature this mass scale must be temperature itself. The quantity u = 
T Jq dr Ai(r, \x\ — > oo) defines the holonomy of the configuration (/? = ^). Due to 
a reflection symmetry u — > — u one can restrict the values of u as < u < 2ttT for 
SU(2). 

Calorons with nontrivial holo nomy are selfd ual configurations that possess BPS 
magnetic monopole constituents | 31 | 32 | 33 | 34 | 35 [ p or an gU(iV) caloron with no net 
magnetic charge there are N constituent monopoles whose magnetic charges add 
up to zero. The masses of the monopoles are determined by the holonomy u and 
thus are oc T. By a recent heroic calculation the one-l oop quantum weight for an 
SU(2) caloron with nontrivial holonomy was derived in^U. 

Since the one-loop effective action of a nontrivial-holonomy caloron scales with 
the spatial volume of the system one is tempted to conclude that these configurations 
do not contribute to the pa rtit ion function of the theory in the thermodynamical 
limit and thus are irrelevant ^Sl. This conclusion, however, is not valid since one can 
show that nontrivial-holonomy calorons are instable under quantum correctionJ^H 
Moreover, if, on spatial average, interacting calorons are shown to be described by 
a quantum mechanically and statistically inert, macroscopic adjoint scalar field <p 
and a pure-gauge configuration a^ 9 then the exponent of minus the associated 
effective action S c i can be factored out in the partition function and thus cancels in 
any physical average. This situation holds even if S c i scales with the spatial volume 
of the system and thus is infinite in the thermodynamical limit. On the microscopic 
level, the generation of (instable) nontrivial holonomy is due to interactions be- 
tween trivial-holonomy calorons mediated by long-wavelength fields that reside in 
the topologically trivial sector of the theory. 

For the SU(2) case it was shown in that the one-loop fluctuations around 
calorons with a small holonomy generate an attractive potential between the two 
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BPS monopole constituents. This implies that monopole and antimonopole anni- 
hilate shortly after they have been created. Thus the likelihood of such a process 
roughly is determined by the finite one-loop quantum weight of a trivial-holonomy 
caloron. Up to an additive correction, which depends on temperature and the 
caloron radius p and which is finite, the effective action equals the classical ac- 

Q 2 

tion S = =%r. For g not too small the likelihood of generating a small holonomy 
is sizable. In the opposite case of a large holonomy a repulsive potential between 
the monopole constituents arises due to one-loop fluctuation Thus monopole 
and antimonopole separate back-to-back, and the caloron dissociates into a pair 
of an isolated monopole and a n isola ted antimonopole which are screened by in- 
termediate caloron fluctuations I 53 | 54 ] B e f ore screening, that is, on the level of the 
classical solu tion the mass of both monopole and antimonopole is much larger than 
temperature ^31. We conclude that the generation and subsequent dissociation of a 
large-holonomy caloron is a very rare process due to an extreme Boltzmann suppres- 
sion. Thus attraction between a monopole and its antimonopole is the dominating 
situation in the ground-state physics of an SU(2) or SU(3) Yang-Mills theory being 
in the electric phase. The macroscopic manifestation of monopole-antimonopole at- 
traction, their subsequent annhihilation, and recreation is a negative ground-state 
pressure. Equating the exponent in the Boltzmann distribution of the monopole- 
antimonopole system before screening with the exponent in the one-loop quan- 
tum weight of the caloron allows for an estimate of the typical distance between a 
monopole and an antimonopole at a given temperature. We will see that, on the 
scale of the inverse temperature, isolated and screened monopoles are very dilute. 
The case of SU(2) has a straight-forward generalization to SU(3): No qualitative 
changes take place in the above discussion when going from SU(2) to SU(3). 

From the selfduality or BPS saturation of the caloron it follows that its energy- 
momentum tensor vanishes identically. Since the macroscopic field <p originates from 
the spatial co rrelations of noninteracting trivial-holonomy calorons of topological 
charge onei^Sl^'s macroscopic energy-momentum tensor vanishes in the absence of 
a coupling to the topologically trivial sector of the theory. This is a derived condi- 
tion which needs to be imposed to fix some of the ambiguities which emerge when 
calculating 0's phase. Namely, one insists on a BPS saturation of the r dependence 
of this phase: A linear second-order equation of motion for </>'s phase can be derived 
from a microscopic definition involving a moduli-space average over a two-point 
correlator, and the requirement of BPS saturation determines the sol utio n up to an 
irrelevant global gauge choice and an irrelevant constant phase shift SSI. 

Subsequently, a potential Ve for the canonically normalized field is derived 
by appealing to the derived information on (jj's phase and to the assumptions that 
a dynamically generated scale exists and that the right-hand side of </>'s BPS 
equation is analytic in (j>. (Here the subscript E refers to the electric phase.) As 
a consequence, cf>'s modulus is oc T~ x / 2 and Ve oc T. Thus the caloron sector 
contributes to thermal quantities in a power-suppressed way at large temperatures. 
Moreover, we will show that the ground state in the electric phase is degenerate with 
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respect to a global electric (SU(2)) and a global electric Z3 (SU(3)) symmetry. 
Therefore, the electric phase is deconfining. 

In the coarse-grained theory a useful decomposition of field configurations 
with trivial topology (only those ones appear as explicit gauge fields) is 



no coupling between Sa^ and a b ^ exists, a fluctuation 5a^ acquires a mass by the 
adjoint Higgs mechanism if [</>, 5a„] 7^ 0. Since the field 4> dynamically breaks the 
gauge symmetries SU(2)-^U(l)and SU(3)— >U(1) 2 , two and six directions in the 
three and eight dimensional adjoint color space acquire mass, respectively. 

We shall discuss why the scale A^, which measures the strength of apparent 
gauge-symmetry breaking by calorons at a given temperature, is physically set at 
a temperature scale Tp where any four-dimensional gauge theory fails to describe 
reality. Common belief is that Tp is comparable to the Planck mass Mp ~ 10 19 GeV. 

How does the existence of the temperature dependent scale \<f>\ influence the 
propagation of gauge modes in the infrared and ultraviolet? Calorons induce quasi- 
particle masses on tree level in the effective theory which are of the order e\<f>\. 
Here e denotes the effective gauge coupling. This coupling measures the interac- 
tion strength between the topologically trivial off-Cartan fluctuations (in unitary 
gauge) and the coarse-grained manifestation of nontrivial topology. Furthermore, 
it is a measure for the screening of the magnetic charge of a monopole. (As far 
as thermodyna mic al quantities are concerned, essentially all excitations are free 
(quasi)particlesSS!.). The evolution of e with temperature follows from the require- 
ment of thermodynamical selfconsistency of this interaction. Except for a small 
range in temperature to the right of the electric-magnetic transition at T c .e, where 



the coupling e is constant: A manifestation of the existence of screene isolated 
and conserved magnetic charges generated by dissociating large-holonomy calorons. 

Infrared cutoffs ~ e\(j)\ arise from a reduction of propagation speed for off-Cartan 
fluctuations by their interactions with calorons, compare with the analogous situa- 
tion for a superconducting material in Fig.[TJ Due to the existence of these cutoffs 
in the loop expansions of thermodynamical quant itie s the probl em of a magnetic 
instability, as encountered in perturbation theory is resolved | 18 | 48 [ j n ^ e m _ 
traviolet, \<j>\ acts as a compositeness constraint by setting a maximal scale for 
the off-shellness of quantum fluctuations. This is a consistent requirement since all 

a The screening of the monopoles is due to surrounding caloron fluctuations of small holonomy and 
not due to Cartan excitations since the latter are not capable of screening static magnetic fields 




e(T) 



log(T - T C , E ) , 



(3) 
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gauge modes, on-shell or off-shell, originate from the nontrivial ground state and 
thus ought not be capable of destroying it. Moreover, plane-wave quantum fluctu- 
ations of an off-shellness larger than \<fi\ 2 are contained, in a coarse-grained form, 
in the pure-gauge configuration a^ 9 . Compositeness constraints are implemented in 
a physical gauge with respect to the unbroken subgroups U(l) or U(l) 2 . The usual 
renormalization program, needed to make sense of ultraviolet divergences in thermal 
perturbation theory, is abandoned in the effective theory: The ground state itself 
provides for a physical ultraviolet cutoff. 

As a consequence of the simultaneous existence of both an ultraviolet and an 
infrared cutoff the contributions of higher loops in the expansion of thermodynam- 
ical quantities are very small. (Tec hnically, an evaluation of two-loop corrections to 
the pressure is quite involved Obviously, the situation outlined so far differs 
substantially from the idea of a Wilsonian flow for nonabelian gauge theories in its 
usual implementation: One derives effective dynamics in dependence of an exter- 
nally set scaleEl k by integrating plane-wave modes harder than k into couplings 
that appear in an ansatz for an effective action. This effective action describes the 
dynamics of gauge modes of maximal hardness k. The dynamics of theses modes 
is, however, not only influenced by integrated-out high-momentum fluctuations but 
also by (spatially) small-scale fluctuations of nontrivial topology. Recall that the 
later can not be expanded in terms of the former because of an essential singularity 
in their weight at a vanishing value of the fundamental gauge coupling. Thus we 
propose an approach which is the converse of the usual picture: Integrate the topo- 
logical sector first and determine subsequently, that is, after spatial coarse-graining, 
what its average effect on the trivial-topology fluctuations is. 

We have already mentioned that this approach to SU(2) or SU(3) Yang-Mills 
thermodynamics implies the existence of three rather than two phases. In the mag- 
netic phase, where the isolated and screened magnetic monopoles of the electric 
phase are massless and thus condensed and where off-Cartan modes are thermo- 
dynamically decoupled, the dual gauge symmetries XJ(1)d or U(l)|, are broken 
dynamically, and the global electric center symmetries Z2 or Z3 are restored in 
the ground state. We will show that the ground state in the magnetic phase is a 
Bose condensate of monopole-antimonopole systems. Each condensate is described 
in terms of a macroscopic and inert complex scalar field and a pure-gauge configu- 
ration. 

A monopole condensate confines fundamentally charged, fermionic and static 
test charges. At the same time, dual and massive gauge modes propagate in the 
magnetic phase. Thus it is appropriate to refer to the magnetic phase as a pre- 
confining phase. The magnetic coupling g, which measures the interaction strength 
between dual gauge modes and condensed magnetic monopoles on the one hand 

b This scale either is continuous, seel2Zlfor a review on gauge theories, or it reflects a scale hierarchy 
originating from the assumed smallness of the coupling constant g at a large temperature T, for 
example k = §tE2] 
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and the screening of center-vortex loops on the other hand, is zero at T c ^e and rises 
rapidly into a logarithmic divergence of the same form as in Eq. ([3]) at a temperature 
T C . M . For SU(2) and SU(3) we have T cM = 0.835 x T C . E and T C . M = 0.877 x T C<E , 
respectively. Thus the magnetic phase occupies only a small region in the phase 
diagram of either theory. This and the fact that the monopole condensates possess 
infinite correlation lengths ~ (monopole mass) - are the reasons why the magnetic 
phase has escaped its direct detection by simulations on finite-size lattices. We will 
show though that it is possible to observe the existence of the magnetic phase in 
a lattice simulations of the infrared insensitive entropy density when using the so- 
called differential method. The latter provides the best-controlled circumvention of 
the infrared problem in simulations on finite-size lattices ^1 

At T Cj m, where g diverges logarithmically, another phase transition takes place. 
All dual gauge modes decouple and thus the monopole condensates, macroscopically 
described by nonffuctuating, BPS saturated complex scalar fields <p (SU(2)) and 
ipx, tf2 (SU(3)) and pure gauges, dominate the thermodynamics. As a consequence, 
the entropy density vanishes at T Cj m and the equation of state is 

p=-P. (4) 

Just like magnetic monopoles are isolated defects in the electric phase there are 
isolated and closed magnetic flux lines in the magnetic phasji which, however, 
collapse as soon as they are created if the magnetic coupling g is finite: Abrikosov- 
Nielsen-Olesen (ANO) vortex loops ^1 Only one unit of flux with respect to U(1)d 
(SU(2)) or either factor in U(l)|> (SU(3)) is carried by a given vortex loop since in 
the electric phase only charge-one calorons dissociate into magnetic monopoles with 
one unit of magnetic chargo This allows for an interpretation of ANO vortex loops 
as center-vortex loops. For SU(iV) the magnetic flux of the latter is determined by 
the differences in phase modulo N of two center elements, see Eq. (TTJ). There is one 
unit of center flux for SU(2), and there are two separate units of center flux for 
SU(3). 

c The fact that only closed loops occur is explained by the absence of isolated magnetic charges in 
a monopole condensate. 

d The core of a vortex line, where U(l)£> (SU(2)) or one factor in U(l)|j (SU(3)) is restored, can 
be pictured as a directed motion of magn etic monopoles (to the right) and antimonopoles (to the 
left) in the rest frame of the heatbatflSSJ see Fig. [2] The magnetic flux, which penetrates a spatial 
hyperplane perpendicular to the direction of monopole or antimonopole motion, is by Stokc's 
theorem measured by a line integral g <f c dz^A® along a circular curve C with infinite radius lying 
in this plane. Here A® denotes the gauge field with respect to the dual gauge group U(l) (SU(2)) 
or either factor in U(l) 2 (SU(3)) generated by the moving chains of monopoles and antimonopoles. 
If we choose to evaluate the line integral in a covariant gauge then the contribution to dz^A® of 
each moving monopole or antimonopole is that of a static monopole or antimonopole since the 
perpendicular part of the gauge field is invariant under boosts along the vortex axis. Thus the state 
of motion of each monopole and antimonopole is irrelevant for its effect on the total magnetic flux 
carried by the vortex as long as the net motion of all monopoles and antimonopoles in a given 
segment defines the direction of the vortex axis: The only thing that determines the magnetic flux 
of the vortex line is the charge of a monopole. 
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Fig. 2. Microscopies of the core of a center vortex. Monopoles (M) and antimonopoles (A) move 
in opposite directions. 

The core-size of a center vortex is determined by the length l g ~ (gl^l) -1 (SU(2)) 
or l g> \ ~ (g\ip\[)~ x , l g< 2 ~ (g\ip2\)~~ 1 (SU(3)) of penetration of the vortex' gauge field 
into a direction perpendicular to the vortex. While center vortices are thick close 
to T c e their core size vanishes at T Cj m- Since the energy of a typical center-vortex 
loorJEH is oc g 1 spin-0 systems composed of a vortex loop and its antivortex loop 
condense at T c ,m- We will show that there is a parameter with discrete values 
characterizing the possible values of a macroscopic, complex scalar field $ which 
describes the vortex condensate. 

The transition to the center phase is of the Hagedorn type and thus nonther- 
mal. An order parameter for this transition is the expectation of the 't Hooft loop 
whose modulus measures the strength of center-vortex condensation. If the 't Hooft 
loop does not vanish then the magnetic Z 2 (SU(2)) or Z 3 (SU(3)) symmetry is 
dynamically broken. These center symmetries are local in four-dimensional space- 
time. Under large U(1)d (SU(2)) or large U(l)|j (SU(3)) gauge transformations the 
macroscopic field $ transforms by multiplications with center elements. 

In the course of the Hagedorn transition the ground state of the magnetic phase 
decays through creation of single and self-intersecting center-vortex loops. In con- 
trast to the magnetic phase, center-vortex loops are stable in the center phase, thus 
are particle-like, and possess a density of states which is over-exponentially rising 
with energy. There are precisely two polarization states for each self-intersecting 
or single loop (spin- 1/2 fermions). After the decay of the monopole condensate is 
completed the new ground state is a Cooper-pair condensate of systems composed 
of a massless single vortex and antivortex loop. The energy-momentum tensor on 
this ground state vanishes identically. This result is protected against radiative 
corrections. 

The center phase is truly confining: A pair of electric, static, oppositely and 
fundamentally charged, and fermionic test charges forms a confining electric flux 
tube because of the presence of condensed electric dipoles (single center-vortex 
loops) in the ground state, and all gauge modes are infinitely heavy. Notice that 
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the absence of propagating gauge modes implies that only contact interactions are 
possible between (thin) center fluxes and the self-intersection points of their vortex 
loops. 

Thermal lattice simulations fail to produce physical results for infrared sensitive 
quantities at temperatures shortly below T c ^e- The center phase as well as the mag- 
netic phase possess infinite correlation lengths. In the center phase this correlation 
length is given by the inverse (vanishing) mass of a single center-vortex loop. Thus 
the fermionic nature of excitations and their over-exponentially rising density of 
states in the center phase escapes simulations performed on finite-size lattices. 

It is self-evident that what was said above has implications for particle physics 
and cosmology. We only would like to mention a few questions that are likely to 
be answered by SU(2) Yang-Mills (thermo) dynamics alone: Electroweak symme- 
try breaking, namely, the origin of the masses of Z and W ± ; the mass hierarchy 
between the two members of a lepton family; the question of whether the neu- 
trino is Dirac or Majorana; the smallness of the cosmological constant on particle 
physics scales; the nature of cosmological and clustering dark matter; cosmic coinci- 
dence; baryon and lepton asymmetries; intergalactic magnetic fields; and large-angle 
anomalies in some of the power spectra of fluctuations in the cosmic microwave 
background. 

The article is organized as follows. In the first part of Sec. [2] we provide prereq- 
uisites on caloron physics. The classical solutions of trivial and nontrivial holonomy 
and their behavior under quantum corrections are discussed for SU(2). As an aside 
we estimate the separation of isolated and screened monopoles in terms of the in- 
verse temperature in the electric phase. The second part of Sec. [5] is devoted to 
the derivation of the macroscopic, adjoint field cj> in terms of a moduli-space and 
S3 integral over the spatial two-point correlations in a caloron-anticaloron system. 
Subsequently, we discuss the full ground-state physics in the electric phase. The 
existence of a dynamically generated scale Ae needs to be assumed to determine 
</>'s modulus at a given temperature. We show that a degeneracy with respect to 
a global electric Z2 (SU(2)) and Z3 (SU(3)) symmetry occurs which proves that 
this phase is deconfining. The last part of Sec. [2] addresses the full thermodynamics 
of the electric phase, including excitations. An evolution equation for the effective 
gauge coupling e is derived and solved numerically. It is observed that an ultraviolet- 
infrared decoupling is manifest in this evolution. We present analytical results for 
the electric screening mass, associated with the massless mode in the SU(2) case, 
and for the two-loop correction to the pressure. 

In Sec. [3] we investigate the magnetic phase. Prerequisites on the BPS monopole 
are given. Considering the average magnetic flux, generated by a screened monopole- 
antimonopole system in a thermal environment, through an S2 with infinite radius, 
where the monopole-antimonopole system is located outside of this S2 , a continuous 
parameter (proportional to the Euclidean time) is derived. This parameter governs 
the temporal winding of a spatially homogeneous, complex, and inert scalar field (p 
in the limit where monopole mass and spatial momentum of the system vanish (Bose 



1, 2008 20:36 WSPC/INSTRUCTION FILE F 



12 Ralf Hofmann 

condensation). For SU(3) two such fields exist, each for every independent monopole 
species. Assuming the existence of a dynamically generated scale Am, the modulus 
of tp is derived. We then discuss the ground-state physics in the magnetic phase. 
It is shown that the electric Z 2 (SU(2)) and electric Z3 (SU(3)) degeneracy of the 
ground state, as observed in the electric phase, gives way to a unique ground state 
in the magnetic phase. Thus we derive test-charge confinement in the magnetic 
phase. Evolution equations for the effective magnetic coupling g are derived and 
solved numerically. A discussion of the full Polyakov-loop expectation (including 
excitations) and an analysis of the critical behavior at the electric-magnetic phase 
boundary are presented. We find that this transition is second order with mean-field 
exponents both for SU(2) and SU(3) with the difference that the peak in the specific 
heat is about three times smaller in the former as compared to the latter case. 

In Sec.[4] we investigate the center phase. We start by providing prerequisites on 
the ANO vortex. In particular, we emphasize the fact that an ANO vortex generates 
negative pressure at finite magnetic coupling g. While a vortex-loop is a particle-like 
excitation at g = 00 it collapses as soon as it is created for g < 00. Collapsing ANO 
or center-vortex loops dominate the (negative) pressure inside the magnetic phase 
where g < 00. In analogy to the monopole condensate we determine a parameter 
(mean center flux through an Si of infinite radius) which governs the expectation 
of a macroscopic, complex scalar field $ describing the Bose condensate of massless 
vortex-antivortex-loop pairs. The values of this parameter are discrete: Two possible 
values for SU(2) and three possible values for SU(3). At T c< m, where g diverges in a 
logarithmic way and where dual gauge modes acquire an infinite mass, the center- 
vortex condensate starts to form under (spin-1/2) particle creation. We construct 
an effective potential Vc for $ involving a scale Ac. We check Vc's uniqueness, 
and discuss how vortex-loop creation takes place by center jumps of $'s phase. An 
estimate for the density of fermion states, created by <&'s relaxation to zero energy 
density and pressure, is provided. As a result, we analytically establish that the 
center-magnetic transition is of the Hagcdorn type. 

In Sec. [5] we discuss how the scales Ae and Am are related by the continuity 
of the pressure across the electric-magnetic phase boundary. We also provide an 
approximate relation between Am and Ac. 

In Sec.[6]we present numerical results for the temperature dependence of thermo- 
dynamical quantities thoughout the electric and the magnetic phase. The following 
quantities are discussed: Pressure, energy density, interaction measure, specific heat 
per volume, and entropy density. While the former quantities are very sensitive to 
the ground-state physics at low temperatures, which is determined by very large 
spatial correlations and thus is inaccessible to finite-size lattices, the entropy den- 
sity is only sensitive to the excitations. This fact makes a comparison of our results 
for the entropy density with those obtained on lattices (employing the differential 
method) useful, all other quantities exhibit quantitative disagreements with their 
lattice-obtained values at low temperatures. 
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In Sec. [7] we discuss implications of our work for particle physics and cosmology. 
We start by addressing the cosmic-coincidence and the old cosmological-constant 
problem in view of a Planck-scale axion, originating from dynamically generated 
and subsequently integrated spin-1/2 fermions at the Planck scale, and in view of 
an SU(2) gauge theory of Yang-Mills scale close to the temperature of the cosmic 
microwave background (CMB) (SU(2) CMB ). This theory is at the electric-magnetic 
phase boundary, and its only massless and unscreened excitation is the photon. 
Throughout cosmological evolution the axion mass is provided by the axial anomaly 
involving nonconfining SU(N) gauge theories of Yang-Mills scales lower than the 
Planck scale. The presence of a Planck-scale axion may explain the particle-number 
asymmetries and CP violation in the weak interactions. Some of the phenomenology 
of the electroweak sector of the SM is addressed in view of leptons being stable 
solitons in the center phase of various SU(2) Yang-Mills theories. These solitons are 
embedded into instable higher-charge excitations with an over-exponentially rising 
density of states which protect their apparent structurelessness seen in scattering 
experiments (the photon couples to the lepton because of mixing) up to center-of- 
mass energies comparable to the mass of the Z boson (with exceptions at momenta 
comparable to the lepton masses). We postdict the mass ratio in terms of 
the mass ratio ^J. Finally, we present some ideas on how fractionally charged 
light quarks and their confinement may arise in Quantum Chromodynamics (QCD) 
in terms of electric-magnetically dual SU(3) gauge dynamics and the fractional 
Quantum Hall effect. 

The last section of the present work briefly summarizes our results. 

2. The electric phase 
2.1. Prerequisites 

2.1.1. The Harrington- Shepard solution (trivial holonoray) 

Calorons of trivial holonomy are the field configurations which enter the definition 
of the phase of the macroscopic adjoint scalar field <f>. We only need to consider the 
SU(2) case since the SU(3) ground-state thermodynamics can be derived from a 
'democratic' embedding of SU(2) calorons. We use the nonperturbative definition 
of the gauge field where the coupling constant is absorbed into the field. 

(Anti) Calorons are (anti)selfdual, that is, their field strength is up to a sign 
equal to their dual field strength 

F lu ,[AVA] = ±F ia ,[AV> A >] (5) 

where the superscript C(A) refers to caloron (anticaloron). Only calorons of topo- 
logical charge one (minus one) enter the definition of 0's phase and thus we will 
focus on this case only. 
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The Harrington-Shepard solutions are given as 

A a 

A%(t, x) = f)a^-j-d u In Il(r, x) or 

A£(r,£) =7? a ^^lnn(r,x) (6) 
where the 't Hooft symbols r\ aiiv and f\ a)iv are defined by 

In Eq. © A a , (a = 1,2,3), denote the Pauli matrices. The periodic solutions in 
Eq. ©, A^' A (0, x) — A^' A (1/T, x) , are generated by a temporal mirror sum of the 
'pre'potential 

n = i + 4 (8) 

X 

of a single BPST (anti)instanton of scale pl^Hlin singular gauge^Sl. Here x 2 = t 2 +x 2 . 
The scalar function n(r, x) in Eq. ^ is given as 

„2 



n(r, a) - ^ 



n— — oo v 



= n(r,r)EEl + ^ 7 ^ / ^ (9) 

£r cosh ^ _ cos (W) 

where r = |x| and (3 = 1/T. Evaluating the integral of the Chern-Simons current 
over a small three-sphere S3, centered at the singular point (r = 0,2? = 0), one 
obtains plus (or minus) one unit of topological charge. For a given value of p the 
solutions in Eq. ([6]) can be generalized by shifting the center from z = to z = (r z , z) 
by the (quasi) translational invariance of the classical action. (The temporal shift t z 
is restricted to < t z < f3 because of periodicity.) In addition, the color orientation 
of each solution can be rotated by global gauge transformations. 

Computing the Polyakov loop at spatial infinity on either of the configurations 
A^ (t, x) and A A (r, x) yields the following result 

r/3 



2 



sinh 



( 

\ J 



P{\x\ — > 00) = V exp 



i I drA^' A (T, \x\ — > 00) 



(i 



1 . (10) 



Thus the Harrington-Shepard solutions possess trivial holonomy. 



2.1.2. The Lee-Lu-Kraan-van Baal solution (nontrivial holonomy) 

For a discussion of (anti)selfdual SU(2) configurations with nontrivial holonomy 
and topological cha rge one (minus one) we use the conventions and closely follow 
the presentation of which to our taste makes the magnetic monopole content 
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z 



Fig. 3. Meaning of the spatial arguments y\,yi entering the solution in Eg. H12I I. The points 
xi,X2 are the core positions of the monopole and the antimonopole. At the point r = 0, a? C m the 
solution is singular. 

most explicit. (In 132 the constituents of nontrivial-holonomy calorons are referred 
to as dyons because the nonabelian magnetic and electric field of each constituent 
is equal and Coulomb-like for large distances away from a given monopole core. 
This property, however, follows from the selfduality of the caloron configuration. 
With respect to the unbroken U(l ) the charge of a constituent monopole is purely 
magnetic (and not dyonic as in^Sj) since the A4 field serves as a Higgs field and not 
as the gauge potential for the electric field.) The existence of these solutions was 
shown by NahmJSl Explicit analytical cons tructio ns were independently performed 
by Lee and Lu and Kraan and van Baal I 32 | 33 [ 

Lee and Lu use antihermitian generators and parametrize the holonomy u as 



where < u < =£-. Using the Nahm data for a monopole coexisting with an an- 
timonopole as an input to the Atiyah-Drinfeld-Hitchin-Manin-Nahm (ADHMN) 
equations (subject to a normalization condition), a selfdual field configuration with 
monopole-antimonopole constituents was constructed in E31 It reads 




(t, |x| -» 00) = -^ A 3 



(11) 



A^x, t) = Ctv^yn u)d + C\v a (y 2 - j - u)C 2 + 



(12) 



where 




(13) 
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Interpreting Va(x \ u) as an adjoint Higgs field, Eqs. (fT3"|) represent the BPS magnetic 
monopole ^1 The matrices C\ , C 2 in Eq. (fT2|) are given as 



Ci = 



Co 



2DN X B{ 
N M 



2DN 2 B 



exp ---s 2 



6XP ( \ ' Sl 



exp I - • s 2 I Q. 



exp -~*2 Q- 



exp ( -*^ rA 3 



TV M 

where Q-t = i(l ± A3) are projection operators. The matrices B 1: B 2 are 



(14) 



B\ = exp 



exp 



B 2 = exp 



exp 



7T 



exp 



A 



IT 



exp 



7T 



exp 



7T 
-I— T 

. 



exp 



2' Sl 



• S2 



A ^ 
2-2 



exp 



exp 



exp 



exp 



A 



S2 



2' S2 



• si 



A ^ 
2' Sl 



and the scalar M. is defined as 

.M = 2 ^cosh si cosh s 2 + yi ■ y 2 sinh si sinh s 2 — cos 
In addition, one defines 



2tt 



TYj = — sinh Si, (i = l,2) 
2/i 



and 



2D 



TV = 1 + — (A^(coshs 2 - (§2)3 sinhs 2 ) + /Y 2 (coshsi + (y 1 ) 3 sinhsi)) 



and 



8 



exp 



■*2 rA3 



(15) 
(16) 
(17) 
(18) 
(19) 



The spatial arguments of the configuration in Eq. (|12p . compare with Fig.[3l are 
defined as 



yi = x-Xi, yi= \yi\ , 
si = uyi , 



s 2 



* = |3| (» = 1,2), 

2tt \ 

-Q-ujm, d = \x 2 -x 1 \. 



(20) 



(A * -sign indicates a unit vector.) Kraan and van Baal show that the distance 
D between the two BPS monopoles can be expressed by the scale p of a trivial- 
holonomy caloron which is deformed to nontrivial holonomy 1221. One has 



(21) 
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Fig. 4. Action density of an SU(2) caloron with nontrivial holonomy plotted on a two-dimensional 
spatial slice. The caloron radius p and therefore the separation D (Eq. 11211 increases from left to 
right while temperature and holonomy are fixed. Figures are taken from a paper by Kraan and 
van Baal. The peaks of the action density coincide with the core positions of the constituent BPS 
monopolcs. 

It was shown in that for y\ <C D 



C 2 ,S~-^= (22) 



and 



D 



A3 cosh 4p — A • yi sinh ^ 



Ci= -. u 2 ===^ +0(1/D). (23) 
^/coshsi - (2/1)3 sinhsi 

Thus C\ is a single-valued unitary matrix, and for 2/1 <C D the configuration in 
Eq. (|12p is an approximate gauge transform of a BPS monopole. Similarily, for 
2/2 D C'2 is a unitary matrix. The difference as compared with C\ for 2/1 <C D is 
that for 2/2 D the matrix C2 induces a large gauge rotation due to an extra factor 



exp 



-ifrA 3 



This gauge transformation inverts the charge of the BPS monopole 
at X2 as compared to the charge of the BPS monopole at x\. There is a singularity 
of the solution at the point (r = 0, a? om ) where 

M -fi + 1 - ^- ^ 2 . (24) 



2tt V 27T, 

This point carries one unit of topological charge. One can show this by expanding 
the solution about (r = 0, a? cm ) and by performing the integral of the Chern-Simons 
current over a small S3 centered at this point. A plot of the action density of a 
nontrivial-holonomy caloron with varying radius p at a fixed temperature and a 
fixed holonomy is presented in Fig.[4j 

On the classical level the masses of the monopoles at x\ and X2 are given as 

mi = 47TM 7712 = 4-7T ( — U J , (25) 

VP / 

respectively. For a large holonomy, that is u ~ 5, we have mi ~ 777,2 ~ 47r 2 T ~ 40 T. 
Thus large holonomy is extremely Boltzmann suppressed. 
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2.1.3. One-loop quantum weights 

In this section we first present the results for the one-loop effective action of a 
trivial-holonomy caloron, which was obtained by Gross, Pisarski, and Yaffe by 
app ealing to the results obtained in I 37 | 38 | 39 l anc j pioneering work of 't Hooft 
Subsequently, we sket ch the results obtained recently by Diakonov, Gromov, 
Petrov, and Slizovskiy SH for the one- loop quantum weight of a caloron with non- 
trivial holonomy. Both results are important for a grasp of the microcopies of the 
ground-state physics in the electric phase. 

Harrington-Shepard solution: 



Th e functional determinant around a Harrington-Shepard caloron was calculated 
in ESI. The result for the quantum weight exp[— S^ B ] of this configuration is given in 
terms of the effective action as 

S eff = ^ + V + WA(a) (26) 

where the dimensionless quantity a is defined as a = 7r4 and 



Ala) = — 
V ' 12 





' d 3 x ({d^) 2 \ 2 f d±x f(d^U ] 



16tt 2 \ IP J J 16tt 2 V II 2 , 



(27) 



The weight exp[— S Bit ] is relevant for the integration over the classical moduli space 
in the presence of one-loop quantum fluctuations. (The nonflat metric is the same 
as in the zero-temperature situation.) 

In Eq. (|27|) the scalar quantities II and LTo are defined in Eqs. ([9]) and ([8]), re- 
spectively. The first integral in Eq. ([27]) is over Si x R 3 while the second integral is 
over R 4 . It is worth mentioning how A(o~) behaves in the high- and low-temperature 
limits a — > oo and a —* 0: 

A(a)^-iloga, (<r-oo), A(a) ^ ~ , (a -f 0) . (28) 

Q 2 

At a given caloron radius p the correction to the classical action thus is large in 
the high-temperature regime, indicating that the contribution of trivial-holonomy 
calorons to the partition function is suppressed, while it is small at low temperatures, 
implying the increasing importance of calorons as the temperature of the system 
drops. The distinction between high and low temperatures is made by a dynami- 
cally generated scale A# which also determines the p dependence of the coupling 
constant g in Eq. J26[) . The latter, by zero-temperature one-loop renormalization- 
group running I 16 | 17 |^ estimatedly becomes larger than unity for p~ x ~ Ae and is 
logarithmically small for p^ 1 ^> A^;. 



Lee-Lu-Kraan-van Baal solution 
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The calculation of the one-loop quantum weight for a caloron of nontrivial holon- 
omy is much harder than for the trivial case. This explains why this result only 
appeared in the literature SU more than six years after the analytical form of the 
nontrivial-holonomy solution was published inGHSH phe expressions are so involved 
that the contribution Z nh of an isolated, quantum-blurred caloron to the total par- 
tition function Z of the theory has so far only been stated in closed analytical form 
in the limit 

D /p N 



p y>1 - (29) 

This, however, is the relevant physical situation, see Sec. 12.1.51 where is is shown 
that (trivial-holonomy) calorons with p^S> f3 dominate the phase of the macroscopic 
adjoint scalar field <f>. As we shall see, the ground-state physics in the electric phase 
is dominated by small-holonomy deformations of the trivial cas e. 

Apart from the restriction in Eq. (|29p the result obtained in SB for Z nh is valid 
for any value of the holonomy, < u < =£■. After the (trivial) integrations of the 

overall color orientation and time translations are performed one obtains S2 

22/3 / p N 5/3 

(2tt + (3uuD)(uD + lj^-^uD + l)^" 1 x 
exp[-V P(u) - 2nDP"(u)} , (D > (3). (30) 



(The number of independent integration variables in Eq. (|30[) is four because 
Jd 3 Xl J d 3 x 2 = 4tt / d 3 x f dDD 2 where x = \{S x + x 2 )). In Eq. ([50]) C ~ 1.0314, 
X\ and Xi are the core positions of the monopoles in the classical solution (compare 
with Fig. [3]), u = =£■ — u, is the Euler constant, V denotes the typical spatial 
volume belonging to the one-caloron system, and A is a scale which is a one-loop 
renormalization group invariant (dimensional transmutation). The functions P(u) 
and P"{u) are given as 



ni \ P 2-2 
P{u) = __ u u . 



p"(u) = 4 



P\ V3 



TT ( 1 

u- - 1 = 



(3\ V3 



(31) 



The function P(u) is always positive for u ^ 0, =^-. The occurrence of the spatial 
volume V in the exponent in Eq. (f30]) would mean total suppression of nontrivial 
holonomy in the naive thermodynamical limit V — > oo. This, however, is not a valid 
conclusion since nontrivial-holonomy calorons are unstable: They either dissociate 
into a pair of BPS monopoles (large holonomy) or they collapse back onto triv- 
ial holonomy by an annihilation of their BPS monopole constituents. This can be 
checked by investigating the second contribution to the exponent in Eq. (|30p . 
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Fig. 5. The quantity -£-P"(u), compare with Eq. H30I I, as a function of the dimensionless holon- 
omy u= ^p. 




Fig. 6. The typical volume spanned by two pairs of BPS monopoles created by the dissociation of 
two calorons whose large holonomy was generated by the interaction of trivial-holonomy calorons. 

In Fig.[5] a plot of -§P"(u) is shown. For < u < |(1 - and for 

+ 773) < u < 2 -| the quantity P"{u) is positive (small holonomy) while it 
is negative in the complementary range (large holonomy) . According to Eq. (|30p 
this means that in the former (latter) case the BPS monopoles experience a linear 
attractive (repulsive) potential. Let us now make an estimate of the typical size of 
an equilateral tetrahedron whose corners are the positions of (screened) magnetic 
monopoles, see Fig.[6l which are generated by the dissociation of a caloron and an 
anticaloron whose large holonomy was created by their interaction. The edge length 
R of the tetrahedron is the typical maximal distance between two BPS monopoles 
generated by a caloron with a holonomy close to maximally nontrivial, w max = 5. 
Once a large holonomy has been created the dissociation of the caloron generates 
the stabilized distance R with probability one. (Once a monopole is at rest with 
respect to the heat bath there is no screening of its magnetic charge by Cartan 
fluctuations but only by small-holonomy calorons in its surroundings. Thus the 
screening of magnetic charge is not described by Eq. (|30p .) Thus it is appropriate to 
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equate the probability for reaching the distance R, where monopoles are sufficiently 
screened to be at rest, governed by Eq. (|30|) , with the thermal probability for excit- 
ing the monopoles in a caloron of large holonomy to start with. Since monopoles 
also are at rest shortly after being created the latter probability is roughly given as 



4tt 2 

exp[-/3(mi + ma)] , mi ~ m 2 ~ — , (32) 

see Eq. (|25[) . Taking only the exponentially sensitive part of the caloron weight into 
account and substituting for V the volume of the tetrahedron, V — jA= R 3 , this 
translates into the following condition: 

7T 2 /R\ 3 2 R o 2 „ 

There exists only a single real and positive solution to this equation. Numerically, 

we obtain R ~ 10.1 /?. So on the scale of the inverse temperature the gas of screened 

magn etic monopoles is dilute. This fits nicely with the lattice results obtained in 
53 54 

While the (extremely small) likelihood for the generation of large-holonomy 
calorons depends on the value of the holonomy only (and not on the distance D) 
this is not true for a caloron with holonomy close to trivial. Since the latter con- 
figuration always collapses back onto trivial holonomy the likelihood for its gener- 
ation is determined by the caloron weight exp[— 5 Gff ] with S Blt given in Eq. (j26|) . A 
strong dependence of S ct{ on D (or p) at a given temperature exists. In contrast to 
exp[— /3(mi + ma)] ~ exp[— 8ir 2 } the weight exp[— S eff ] is sizable at S e{t 's minimum 

"min • 

We conclude that attraction between a BPS monopole and its antimonopole 
(small holonomy), which are in equilibrium with respect of their creation and anni- 
hilation, by far dominates the ground-state physics as compared to the case where 
monopole and antimonopole repulse one another (large holonomy) . Macroscopically, 
this situation expresses itself by a negative pressure of the ground state. We shall 
compute the temperature dependence of this pressure in Sec. 12. 1.61 



2.1.4. Microscopic definition for the phase of an adjoint and macroscopic 
scalar field <fi 

The results that were discussed in the last two subsections are important for an un- 
derstanding of the infrared physics in the electric phase. The detailed microscopic 
dynamics is very complicated and, as it seems, it is impossible to derive macroscopic 
quantities such as the pressure or the energy density or the mass of thermal quasi- 
particles by performing literal ensemble averages on the microscopic level. What 
turns out to be feasible and thermodynamically exhaustive is to compute the spa- 
tial average (spatial coarse-graining) over the physics generated by the topologically 
nontrivial sector. This procedure introduces the concept of a macroscopic, thermal 
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ground state. As far as thermodynamics is concerned one still obtains exact results 
this way. The advantage of such an approach is that the complications of a mi- 
croscopic calculation are avoided. Once the ground-state physics is understood and 
quantitatively described its effect on the propagation of trivial-topology modes can 
be investigated. 

If the ground state is to be characterized by a macroscopic field other than 
a pure-gauge configuration then, by spatial isotropy, this macroscopic field must 
be a Lorentz scalar <f>. Moreover, in a pure Yang-Mills theory, where all local fields 
transform under the adjoint representation of the gauge group, the composite field <p 
needs to transform in an adjoint way under the remnants of a microscopic, spacetime 
dependent gauge transformations. Since space dependent gauge transformations are 
constant on the macroscopic level (due to the spatial average) no space dependence 
of 4> occurs in any chosen gauge. Apart from its modulus, which is governed by a 
dynamically emerging scale A^ and temperature, the only nontrivial information 
on 4> is the r dependence of its color orientation in a given gauge. In the following 
we will refer to </>'s color orientation as ^>'s phase. 

Let us imagine a (hypothetical) Yang-Mills world where the only field configura- 
tions allowed to contribute to the partition function are classical and noninteracting 
caloron configurations of trivial holonomy. We will show in Sec l2.1.5l that it is con- 
sistent to adopt this point of view in the derivation of the macroscopic ground-state 
physics. Since the Yang-Mills scale Ae can not be computed we focus on the com- 
putation of </>'s phase first. Because </>'s phase is a ratio of the field and its modulus 
and hence dimensionless the associated me asure for the p-average is flat. 

We closely follow the presentation in u for the remainder of this section and for 
Sec. l2.1.5l Due to the selfduality of calorons any local definition of <f>'s phase yields 
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the trivial result zero. Thus we start by defining: 



\4>\ 



(r) ~ tr 



-F^[A a (p,f3)}((T,0)) {(T,0),(T,x)}[A a ( P ,(3)} x 



/3°1! / d 3 x I dp 

T 

F„ u [A a (p, P)] (( T , x)) {(t, x), (r, 0)} [A a (p, /?)] + 
-1 2! /" d 3 x /" d 3 y y dp 

yF M [4(p^)]((r,0)) {(T,0),(r,f)}[A Q (p,/3)] x 
F Xv [A a (p,p)]((T,x)) {(T,x),( T ,y)}[A a (p,f3)] x 
F„ M [A Q (p,/?)] ((r,y)){(r,y),(r,0)} [A Q (p,/3)] + 

/T 2 3! y d 3 x J d 3 y J d 3 u J dp 

y^[4„(^)]((r,0)) {(T,0),(T,x)}[A a (p,/3)] x 

F Xv [A a (p,P)]((T,x)) {(T,x),(T,$)}[A a (p,f3)] x 

F„„[i4«(p,0)] ((r,y)){(T,y),(r,M)} [A a (p, 0)}F Kfl [A a (p, 0)] ((r,u)) x 



{(r,u),(T,0)}[ J 4 Q (p, / 9)] + 



(34) 



A number of comments are in order: The dots in (|3"4")) stand for the contributions of 
higher n-point functions and for reducible, that is, factorizable contributions with 
respect to the spatial integrations. The factors (n — 1)! in front of the n-point con- 
tribution measures the multiplicity of the corresponding integral. Factors (3 n ~ 2 are 
needed to make the contribution dimensionless. The argument A a (p, (3) (spacetime 
dependence suppressed) refers to either a caloron or an anticaloron configuration, 
the Harrington-Shepard solutions of Sec. 12.1. II Moreover, the following definitions 
apply: 



|0| ^ - tr0 2 , 
{(T,0),(r,x)}[A a ] =V exp 

{(r,f),(r,0)}K] =7>exp 



(r,x) 

' / dy p Ap{y,p) 

(t,0) 
(t,x) 

-' I dy(3Ap{y,p) 

(t,0) 



(35) 



where V is the path-ordering symbol. Under a microscopic gauge transformation 
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Q(y) we have: 

{(r, 0), (r, x)} [A a ] -» fit (( T , 0)) {(r, 0), (r, f)} [A a ] Q((r, x)) , 
{(r, g), (r, 0)} [AJ - nt(( Tj f) ) { ( Tj £), ( T) o)} [A a ] n((r, 0)) , 
F^[A a ]((r,f))^nt((T,x))^K]((r,f))fi((r,f)), 
F^L4 Q ] ((t, 0)) - nt(( Tj o)) i^L4 a ]((r, 0)) fi((r, 0)) . (36) 

As a consequence of Eq. (f36|) the right-hand side of ([34| transforms as 

W\ {T) ^ Rab{T) W\ {T) (37) 

where the SO (3) matrix R a b{r) is defined as 

R ab (T)\ b = n((T,o))\ a n\(T,o)). (38) 

Thus we have defined an adjointly transforming scalar in (|34j) . In addition, we have 
just shown that only the time-dependent part of a microscopic gauge transformation 
survives on the macroscopic level. (Shifting the spatial part of the argument (r, 0) — > 
(t,z) in (f34|) introduces a finite parameter z to the gauge rotation i? ah : R ab (r) — > 
R ab (T, z). Such a shift, however, introduces an arbitrary but finite mass scale 
into the definition of 0's phase which, on the classical level, is absent. Also, a finite 
value of \z\ would introduce an explicit breaking or rotational symmetry into the 
definition (l34|) . Thus we have z = 0. Moreover, the integration path connecting 
the points (r, 0) with (r, x) in Eq. (I35|) ought to be a straight line since a spatial 
curvature would imply the existence of a mass scale other than temperature on 
the classical level.) Integrations over shifts r — > r + r s (0 < r s < /?) project a 
nontrivial (periodic) r dependence of 0's phase onto zero and thus are forbidden. 
Integrations over global gauge rotations are forbidden for the same reason. Spatial 
shifts x — > x + x s ,y — > y + x s , - ■ ■ leave the integrals in (|3"1)) invariant. These 
averages are already performed. Thus the only admissible integration over moduli- 
space parameters is over p with a flat measure. 

In (|34|) the ~ sign indicates that both left- and right-hand sides satisfy the same 
homogeneous evolution equation in r 

v [w\\ =0 - m 

Here T> is a differential operator such that Eq. (|39|) represents a homogeneous dif- 
ferential equation. As it will turn out, Eq. (|39p is a linear second-order equation 
which, up to global gauge rotations and a choice of winding sense, determines the 
first-order or BPS equation whose solution 4>'s phase is. (The ambiguities in the 
evaluation of the right-hand side span the solution space of Eq. (|39|) . and thus V is 
uniquely determined by (|34jl .) 

We now discuss why n-point functions with n > 2 do not contribute to the right- 
hand side of (|34|) . Since the classical (anti)caloron action S = t- and the classical 
moduli-space metric are independent of temperature we conclude that no explicit 



1, 2008 20:36 WSPC/INSTRUCTION FILE F 



Nonperturbative approach to Yang-Mills thermodynamics 25 

dependence on may occur in the definition of <f>'s phase. For n > 2, however, 
explicit (3 dependences do occur, see (l34l) . We conclude that these contributions to 
<^'s phase do not exist. 

What about calorons of higher topological charge? Some of these solutions have 
been constructed, see for example | 43 | 44 | essential difference to the charge- 
one case is that more dimensionful moduli occur than just the parameter p. For 
example, for charge-two configurations there is a spatial core separation between 
the two seed instantons and an additional instanton radius p' . The reader may now 
convince himself that along the lines of (fM)) a nonlocal definition of </>'s dimensionless 
phase, which would also have to include integrations over the additional moduli of 
dimension length, is impossible for higher-charge calorons. (This is certainly true for 
the integral over two-point functions. For every increment in n there is an increase 
in power of length scale by one unit arising from an additional d 3 x x F^ v . This 
makes the situation even worse in comparison to the two-point case.) 

2.1.5. Essentials of the calculation 

Before we dive into the essential parts of the calculation, which will lead to the 
unique determination of the operator V in Eq. (|3"§|) , we would like to discuss a 
condition which severely constrains the possible solutions to this equation. 

By (anti)selfduality the energy-momentum tensor vanishes identically on a 
caloron or an anticaloron, 

9^[A^}=0. (40) 

Since </f>'s phase is obtained by an average over (the admissible part of) the moduli 
space of a caloron-anticaloron system (no interactions) the macroscopic energy- 
momentum tensor # M „[</>] should vanish identically as well, 

0uM = 0- (41) 

In a thermal equilibrium situation, described by Euclidean dynamics, this is true if 
and only if the r dependence of 4> (or <j>'s phase) is BPS saturated. Thus cj) solves 
the first-order equation 

= v£ /2) (42) 

where Vg denotes the 'square-root' of a suitable potential. (The fact that 
an ordinary and not a covariant derivative appears in Eq. (|42|) is, of course, 
tied to our specific gauge choice. If we were to leave the (singular) gauge for 
the seed (anti)instanton, in which the solutions of Eq. © are constructed, by a 
time-dependent gauge rotation O(r) then a pure-gauge configuration A p ^ 9 '(t) = 
i8 ^Cft d T Cl would appear in a covariant derivative on the left-hand side of Eq. (|42[) . 
Also recall the fact that the heat bath breaks boost invariance. This is encoded in 
the noninvariance of Eq. (@2j) under 0(4) rotations.) V E = tr (V^ 1/2) ) V^ 1/2) . As 
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we will see below, the right-hand side of Eq. (|42p is determined only up to a global 
gauge rotation and a choice of winding sense. 

Let us now discuss essential details of the calculation of the right-hand side of 
(|34| which, after what was said in Sec. 12. 1.41 reduces to 



W\ (T) 



d 3 x 



dp^-F^[A a (p,/3)}((T,0)) {(T,0),(T,g)}[A a (p,f3)] x 



tr 

F^[A a (p, 0)] ((r, x)) {(r, x ), (r, 0)} [A a (p, 0)] 



(43) 



where a sum over the contributions of a trivial-holonomy caloron and an anticalorons 
is to be performed. 

Since the integrand in the exponent of the Wilson line {(r, 0), (r, x)} [A^' A ] is a 
hedgehog the path-ordering prescription can be omitted. For the caloron contribu- 
tion one obtains 



dp 



crx — x 
r 



(d 4 n(T + T C ,o)) 2 

n 2 (r + r c ,0) 
<9 r <9 4 n(T + t c , r) 



2 djn(T 

3 n(r- 



T C ,0) 



|4cos(2.g(r + t c , r)) x 



rc,0) 

(a r n(r + r c ,r))(a 4 n(r + T C ,r)) 



II(r + r c ,r) 
sin(2fl-(r + r c ,r)) 



II 2 (r + r c ,r) 
(5 4 n(r + T C , r)) 2 - (a r n(r + r c , r)) 2 



I1 2 (t + tc, r) 



5 2 n(r + r c , r) - c) 2 n(T + T C , r) 
LT(t + r c ,r) 



where 



l(r + Tc,r) = 



ds — di In I1(t + rc, sr) , 



(44) 



(45) 



the function LT(r, r) is defined in Eq. (|9|), and tq refers to a constant but arbitrary 
temporal shift of the caloron center (0 < rc < (3). It is worth mentioning that the 
integrand in Eq. (|45|1 is proportional to S(s) for r f3. The dependences on p and 
/3 are suppressed in the integrands of (|44|) and Eq. ([45|) . 

As compared to the contribution of the caloron there are ambiguities in the 
contribution of the anticaloron: First, the r dependence of the contribution of the 
anticaloron may be shifted by ta (0 < ta < /3). Second, the color orientation of 
caloron and anticaloron contributions may be different. Third, the normalization 
of the two contributions may be different. To see that this is true, we need to 
investigate the convergence properties of the radial integration in (|4*4"| . It is easily 
checked that all terms give rise to a converging r integration except for the following 
one: 

9 2 n(r + r c ,r) 



2— sin(2 5 (r + rc,r)) 

r HXt + tc,?") 



(46) 
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x 2 




Fig. 7. The axis for the regularized azimuthal integration. 



Namely, for r > R ^> (3 (l46)l goes over in 



x a np 2 sin(2g(r + rc,r)) 
r /3r 3 



(47) 



Thus the r-integral of the term in (|46[) is logarithmically divergent in the infrared: 
(The integral converges for r — > 0.) 



Recall, that g(r + Tc,r) behaves like a constant in r for r > R. The angular 
integration, on the other hand, would yield zero if the radial integration was regular. 
Thus a logarithmic divergence can be cancelled by the angular integral to yield some 
finite and real answer. To investigate this in more detail, both angular and radial 
integration need to regularized. 

One introduces a regularization, conveniently we have chosen dimensional reg- 
ularization in with a dimensionless regularization parameter r)c > 0, for the 
r-integral in Eq. (|48[l while the angular integration can be regularized by introduc- 
ing defect (or surplus) angles 2r/ c in the 9 integration (azimuthal angle in the X1X2 
plane) . Any other plane for the azimuthal angular integration could have been cho- 
sen. Moreover, the value of ac is determined by a (physically irrelevant) initial 
condition, as we will show below, see Fig. [71 Together, the choice of the regulariza- 
tion plane and of the angle ac amount to a global choice of gauge: an apparent 
breaking of rotational symmetry by the angular regularization is nothing but a 
gauge choice.) Without restriction of generality (global gauge choice) we may also 
for the contribution of the anticaloron use an axis for the angular regularization 




(48) 
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which lies in the x 1X2 plane, but with a different angle a a- Then we have 



\<t>\ 



±^c (S a i cos a c + S a 2 sin a c ) A 



2tt(t + tc) 



2tt(t + t a ) 



±Sa (5 a i cos a a + Sa2 sin ola)A\ 

\ P J 

(49) 

where Hp , Ha, tc, ta, ac, and a a are undetermined. (He, Ha are the ratios of 
rjc,A and f7 Cj4 , respectively.) The function A ( ^ is defined as 



.4 



2ttt 
T 



32 7T 7 

y /? 



I d P Jim sin(2g(r,r)) 



ttV + /3 2 (2 + cos(^ 



2vrV + /? 2 (l-cos(^)) 



(50) 



Eq. (|49|) and Eq. (|50|) provide the basis for fixing the operator T> in Eq. (f39|) . To 
evaluate the function .A ( 2g£ ) in Eq. ([50)) numerically, we introduce the same cutoff 



for the p integration in the caloron and anticaloron case as follows: 

dp^ J dp, (C>0). (51) 

This introduces an additional dependence of A on (. In Fig. [S] the t dependence of 
A for various values of £ is depicted. Therefore we have 

j- ~ 272 C 3 ( H c (<5 Q i cos a c + S a 2$inac) sin \^-(t + t c 



\'P\ 



(3 



-E A (S a i cos OA + S a2 sina^) sin ( — (r + t a ) 



2vr 



= 4> a ■ (52) 

Just like the numbers He and Ha are undetermined on the classical level due to the 
invariance of the classical action under spatial scale transformations so is the number 
£. It is clear, however, from Eq. ([52]) that the integral in A is strongly dominated by 
p- values close to the upper integration limit. Let us now discuss the physical content 
of (|52|) . For fixed values of the parameters ( 3 He, C 3 Ha, ^ and ^ the right-hand 
side of Eq. (|52[) resembles an elliptic polarization in the X1X2 plane of adjoint color 
space. For a given polarization plane the two independent numbers (normalization 
and phase-shift) for each of the two oscillations parametrize the solution space of 
the second-order linear differential equation 



V4> = 0. 



(53) 
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^Jl/|3) T 



Fig. 8. A as a function of -^-t for £ = 1,2,10. For each case the dashed line is a plot of 

max .A X sin (^f" T )- We have fitted the asymptotic dependence on f of the amplitude of A as 

A {^jt = f ,C) = 272£ 3 , (C > 10). The fit is stable under variations of the fitting interval. For 
the case £ = 10 the difference between the two curves can not be resolved anymore. 



From ([52]) we observe that the operator T> is 

'2vr 



V = di 







(54) 



The ambiguities in Eq. (|52"|) parameterize the solution space of Eq. ([53]) for a given 
polarization plane which depends on a global choice of gauge. Thus the differential 
operator T> is uniquely determined by Eq. (|52"|) . What is needed to assure the valid- 
ity of Eq. (|4Tj) is a BPS saturation of the solution to the linear Eq. (|53|) since the 
modulus of <fi may not depend on r in thermal equilibrium. 

Thus we need to find first-order equations whose solutions are traceless, hermi- 
tian and solve the second-order equation ([53} . The relevant two first-order equations 
are 

2iri 



d r4> = ± — A 3 </>. 



(55) 



Obviously, the right-hand sides of Eqs. (|5"5"|) are subject to a global gauge ambiguity 
(associated with the choice of polarization plane in which the regularization of the 
azimuthal angular integration is carried out) and a choice of sign: Any normalized 
generator other than ±A 3 could have appeared. Moreover, a solution to either of 
the two equations ([53)) also solves Eq. (p)3"|) for a given polarization plane, 



7T 



2ni 2ni 

T T 



Solutions to Eqs. ([55]) are given as 
4> — C Ai exp 



27ri 
T— A 3 (r- 



ro) 



(56) 



(57) 
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where C and To denote real integration constants which both are undetermined. We 
set To = in what follows. The solutions in Eq. (|57[) represent a circular polarization 
in the x\X2 plane of adjoint color space and thus indicate that the field <j> winds 
along an S\ on the group manifold S3 of SU(2). Both winding senses appear but 
can not be distinguished physically: A change in winding sense does not affect the 
potential nor does it affect the admissibility of the transformation to unitary gauge, 
see Sec. l2~2T2"l 

2.1.6. <p's modulus and potential 
SU(2) case: 

The information in Eq. (|5T[) on <^'s phase can be used to infer its modulus once the 
existence of an externally given mass scale is assumed. (The scale A^ determines 
the typical distance between caloron centers at a given temperature.) As long as 
no interactions between trivial-holonomy calorons are allowed for this is consistent 
since the BPS saturation of <f> forbids the occurrence of (gravitationally) measurable 
effects: The macroscopic energy-momentum tensor vanishes identically, and thus 
assuming the existence of the scale A^ docs not yet influence the ground-state 
physics. We have 



In order to reproduce the phase in Eq. (|57| a linear dependence on (f> must appear 
on the right-hand side of the BPS equation (|42p . Moreover, this right-hand side 
ought not depend on (3 explicitly and must be analytic in <f>. The former require- 
ment derives from the fact that <j) and its potential V are obtained by functionally 
integrating over the (admissible part of the) moduli space of a caloron-anticaloron 
system with no interactions. The associated part of the partition function does not 
exhibit an explicit (3 dependence since the action and thus the weight are (3 inde- 
pendent on the moduli space. Thus a (3 dependence of V or V^ 1 / 2 ) can only be 
generated via the periodicity of <f> itself. The latter requirement derives from the 
demand that the thermodynamics at temperature T + ST to any given accuracy 
must be smoothly derivable from the thermodynamics at temperature T for ST suf- 
ficiently small provided no phase transition occurs at T. This is done by expanding 
the right-hand side of the BPS equation (finite radius of convergence) which, in 
turn, is the starting point for a perturbative treatment with expansion parameter 

ST 



Linearity, analyticity, and no explicit dependence of (3 only allow the BPS equa- 
tion for <j> to be one the two following possibilities: 




(58) 



d T 4> = ±i A E A 3 <j> 



(59) 



or 



d T <f> = ±iA% A 3 (T 1 



(60) 
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where 4> = rfa- Recall that 



n=0 



has a finite radius of convergence. According to Eqs. (|58[) and (|57p we may write 

/ 27T? \ 

- |0|(/3,A B ) x Ai exp f T— A 3 tJ . (62) 



Substituting Eq. J62J) into Eq. ([59]) yields 



A £ = ^ (63) 



which is unacceptable since is a constant scale. For the other possibility Eq. (|60p . 
we obtain 



pA| . / A| 



when substituting Eq. (f6"2"]) into Eq. ([BH]) . This is acceptable and indicates that at 
T ^S> Ae <f>'s modulus is small. The right-hand side of Eq. (ffJOf defines the 'square- 
root' F (1/2) of a potential V(|<£|) = tr (V^ 2 ^ V {1 W = A|tr0- 2 . The equation 
of motion (|60|) can be derived from the following action: 

5 = tr J dr J d 3 x (d T 4>d T (l) + k E 4>~ 2 ) . (65) 

Notice that due to BPS saturation it is not possible to add a constant to the 
potential in Eq. (|65|) without changing the ground-state physics. (In fact, adding a 
constant, the modified BPS equation would not admit periodic solutions anymore.) 
The scale \<fi\ must be interpreted as the maximal resolution that remains after 
the spatial coarse-graining over calorons and anticalorons is performed. As we shall 
show later, a critical temperature 2itT c .e = 13.867 A^ exists. Thus, expressing the 
critical cutoff l^l -1 = J A:i 2 ^ e in units of c> e, yields 8.22; for T > T c ^e this 

number grows as (T/T Ci e) 3 ' 2 - But cutting off the p- and r-integration at > 8.22/3 
perfectly represents the infinite- volume limit in Eq. (|50[) ! 

The ratios of the mass-squared of 0-field fluctuations, dPn V(|0|), and the com- 
positeness scale \<f>\ squared or T 2 are given as 

-gj5- = 12A|, J<±-=48**, (66) 

where \e = We will show in Sec. 12.3.41 that A^ > 13.867 in the electric phase. 
Thus both ratios in Eq. (|66| are much larger than unity: The field <j> is quantum 
mechanically and statistically inert. It represents a background for the dynamics of 
the topologically trivial sector after spatial coarse-graining. As a consequence, our 
assumption that only noninteracting calorons of trivial holonomy contribute to the 
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average in Eq. (|13"]) is consistent. 
SU(3) case: 

For SU(3) we write three sets of SU(2) generators as 

/010\ /0-i0\ 

Ai = [ 1 I , A 2 = ( i I , A 3 = | -10 | , (()?) 

and 

Ai = [ | , A 2 = | ) , A 3 = | | , (68) 

and 

Ai = I 1 I , A 2 = I -i I , A 3 = I 1 I . (69) 












One generator is dependent. This just refle cts the fact that the group manifold of 
SU(3) locally is not S3 x S3 x #3 but S 3 x S 5 1 59 | 60 |_ ^ se ^ f independent generators 
is obtained by replacing the two matrices A3 and A 3 by the single matrix 

A 8 = 4=(a 3 + A 3 ) = 4t I 1 1 (70) 




V3 V M V V3 

and by keeping the other matrices. The result is the familiar set of Gell-Mann 
matrices generating the group SU(3). 

For the case of SU(3) the field <f> may wind in each of the above SU(2) algebras. 
Except for the points t = 0, ^, where it jumps into a new algebra, a solution to 
the BPS equation 

A 3 ^, (0<r<|) 



A 3 ^, (l< <) 



^ = ±^1^3^, (|<T<f) (71) 



is given as 



Ai exp (t^A 3 t) , (0 < r < 



*k \ Axexp^AaCr-DK (f < r < f ) ,72, 
A ie xp(T^Aa(r-f)) , (f < r < 0) . 

A 6 

Notice that the potential Vg = 2y^- is the same on the configuration 4>{t) in 

Eq. (|72[) as for the SU(2) case and that by the same calculation one shows its 
quantum mechanical and statistical inertness. 
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2.2. A macroscopic ground state 

The action Eq. (|65|) governs the dynamics of 4>- We have not yet included caloron in- 
teractions, mediated by the topologically trivial sector, which change the holonomy 
of calorons and induce interactions between their (BPS monopole) constituents. 
This is the objective of the present section. 

2.2.1. Pure- gauge configuration 

The action Eq. (|6"5")) can be extended to include topologically trivial configurations 
0^. This is accomplished by a minimal coupling d T (f> — > d^<p + * e [0, %i] = D^ifi and 
by adding a kinetic term for these configurations. Here e denotes the effective gauge 
coupling. The total Yang-Mills action S, governing the electric phase, can thus be 
rewritten as 

S = ix J drj d 3 x {^G^Gau + D^D^ + A%<j>-^j , (73) 

where G uv = G« v % and G« v = dy,a% - d v a^ + e f ahc a h ^a c v . The classical equation 
of motion for a M , derived from the action (|73p . reads 

DpG^ = ie[4>, D v <\>\ (74) 

There is a pure-gauge solution ajf to Eq. (f?4"|) with D v <f> = 0. Thus the total action 
density of the ground state ((f>,a b l f) is the potential Ve = A.itA E T, corresponding 
to an energy-momentum tensor M „ = Ve 5^ or P gs = —p 9S = — 4tt A e T (P gs , p 9S 
the ground-state pressure and energy density, respectively): The so-far hidden scale 
Ae becomes (gravitationally) measurable by coarse-grained interactions between 
calorons. 

SU(2) case: 

In the background 

m = / J Al exp (t^A 3 t) (75) 

we have 

= ±6^ A 3 . (76) 

SU(3) case: 

In the background qb of Eq. (|72p the pure-gauge solution to Eq. (|74[) with D„qb — 
reads 

r a 3 , (o<r<§) 

< 3 = ±^-A As, (|<r<f) (77) 
Gfj ( As, (f <r</3). 
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2.2.2. Polyakov loop and rotation to unitary gauge 

Here we would like to investigate whether the ground state, described by (0, a^ 9 ), is 
degenerate with respect to the global electric Z 2 (SU(2)) or Z 3 (SU(3)) symmetry 
under which the Polyakov loop P transforms as P — > ZP where Z £ Z2 (SU(2)) 
or Z e Z3 (SU(3)). We will refer to the gauge, where 4>'s phase is r dependent as 
in Eq. ([62"]) or in Eq. (j72j) . as winding gauge. (Microscopically, this is the singular 
gauge for an instanton in which the Harrington-Shepard solution is constructed.) 



SU(2) case: 

Evaluating the Polyakov loop on the configuration a* 9 of Eq. ([76]) , we have 

P[aJ»] = exp [±in\ 3 } = -1 2 . (78) 

We are searching a gauge transformation f2 €SU(2) to the unitary gauge <j> = 
|0| A3 and a^ 9 = 0. A periodic but not differentiable gauge transformation doing 
this can be obtained from a nonperiodic but smooth gauge transformation fi by 
multiplication with a local center transformation Z and by multiplication with a 
global gauge transformation £l g f. 

q = n(T)z( T )n gl , (79) 

where fi(r) = exp[T«7r^A 3 ], Z(t) = (28(t - f ) - l) 1 2 , and fl g i = expHf A 2 ]. 9 
denotes the Heavyside step function: 



0, 


(x 


<0), 


1 

2 • 


(x 


-0), 


1. 


(x 


>0). 



e(x) = i i (a = 0) , . (80) 



Applying fi to = a^ 9 + <5a M , where <5a M is a periodic fluctuation in winding gauge, 
we have 

= (nrg^^V (si) 

Since QQ g i(0) — ~tttt g i(P) we conclude that if the fluctuation 6a^ is periodic in 
winding gauge then it is also periodic in unitary gauge. It thus is irrelevant whether 
we integrate out the fluctuations 5a u in winding or unitary gauge in a loop expansion 
of thermodynamical quantities: Hosotani's mechanism 1 ^- 1 does not take place. What 
changes though under the gauge transformation (l is the Polyakov loop evaluated 
on the background configuration a* 9 : 

P[4 9 ] = -1 2 ^P[< 9 =0] = 1 2 . (82) 
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We conclude that the theory has two equivalent ground states and that the global 
electric Zi symmetry is dynamically broken. Thus we have shown that the elecric 
phase is deconEning. 

SU(3) case: 

Let us now compute the Polyakov loop on the configuration a^ s of Eq. ([77)) . We 
have 



exp 
exp 



±z-A 3 exp 
^|(±A 3 ±A 3 ±A 3 ) 



exp 



(83) 



The + or — sign can be chosen independently for each SU(2) algebra. The following 
combinations are possible: 

± (+A 3 + A 3 + X 3 ) = ±2 A 3 , 



± -A 3 - A 3 + A 3 = ±2 A 



± 



± 



(-A 3 + A 3 + A 3 ) 



±2 A, 



-As - A 3 + A 3 = 



(84) 



The corresponding values of the Polyakov loop are 



P ± 




P$ = ( Oexpfi^ 1 




(85) 



P4 is a trivial representation of the center group. The set Pf, P^, P3 closes under 
multiplication with the center elements P = exp[±^p] 1 3 , 1 3 . It is a six dimensional, 
reducible representation of the center group. The two three dimensional irreducible 
representations, which collapse on one another, are spanned by 
1 



rla Pf 



1 , 2wi, 
gexp[T^-]l3 



]1 3 Pf 



P^ 



(86) 



We conclude that the ground state has a Z 3 degeneracy: The electric Z3 symmetry 
is dynamically broken and thus we have discussed a deconfining phase. 

What about a gauge rotation to unitary gauge a* 9 = and <fi = | <^> | A3 or <f) = 
|0|A 3 or <p = |^| A3? Such a gauge transformation Q is given as 



exp[=Rf t\ 3 ] exp [-if A 2 ] , 
<> = ( exp^f (r - §)A 3 ] exp [-if A 2 ] , 
exp[T»f(r- f )A 3 ]expHfA 2 ], 



(0<r<§) 
(§<r<f) 
(f<r<0). 



(87) 



1, 2008 20:36 WSPC/INSTRUCTION FILE F 



36 Half Hofmann 

By construction f2 is periodic, at r = (3 it jumps back to its value at r = 0, and 
thus a fluctuation <5a M , which is periodic in winding gauge, is also periodic in unitary 
gauge. 

A comment concerning SU(iV) theories with TV > 4 is in order. We only discuss 
the case when TV is even. Since at any r the maximal number of independent SU(2) 
subgroups contributing with calorons of topological charge one to the macroscopic 
field cj> is N/2 the dynamical gauge-symmetry breaking is not as maximal as it is 
for SU(2) and SU(3). For example, SU(4) breaks only down to SU(2) 2 xU(l). The 
question is whether there is a single critical temperature T Cj e where the unbroken 
nonabelian subgroups get broken to U(l) factors by the condensation of color mag- 
netic monopoles into macroscopic adjoint Higgs fields and pure gauges, and where 
the abelian magnetic monopoles with respect to the remaining U(l) factors con- 
dense as soon as they are created, or whether this is a stepwise process. The lattice 

IOOI 

seems to favor the former situation 12 ^ which, however, does not appear natural to 
the author. 

2.3. Excitations 

Now that the derivation of a macroscopic ground state for the electric phase is com- 
pleted we are in a position to discuss the properties of its on-shell excitations and the 
role of residual quantum fluctuations. The temperature dependent mass is a mea- 
sure for the strength of coupling between gauge modes and the nontrivial ground 
state. The evolution of this coupling with temperature is a manifestation of the 
thermodynamical selfconsistency of the separation into topological configurations 
and gauge modes with trivial topology. As we shall see, this evolution represents a 
decoupling between high and low temperature physics and, for temperatures suffi- 
ciently above the critical temperature T Cj e, indicates the conservation of magnetic 
charge associated with the isolated and screened BPS monopoles that are liberated 
by the dissociation of calorons with a large holonomy. 

2.3.1. Mass spectrum of thermal quasiparticles 

We first discuss some general aspects of the mass spectrum and then specialize to 
the cases SU(2) and SU(3). We refer to gauge modes which acquire a quasiparticle 
mass on tree level by the adjoint Higgs mechanism as tree-level heavy (TLH) and 
to those which remain massless as tree-level massless (TLM). 
In unitary gauge the mass spectrum calculates as 



777' 



2 



2e 2 tr [cj),t a }[<f>,t a 



(88) 



a 



where t a are the group generators normalized as trt a t b 
diagonal generators can be chosen as 



\5 ab . For SU(iV) off- 



t: 



H = 1/2 (Stf + Stf) , t£ = -i/2 {SlSj 
{I = !,■■■ ,N; J > I; r,s = !,■■■ ,N) . 




(89) 
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This yields 

mjj = mjj = e 2 (</>,- 4>jf (90) 
where 4>i, 4>j denote the diagonal elements of <fi i n unitary gauge. 
SU(2) case: 

In this case <f> breaks the gauge symmetry dynamically down to U(l). Thus we have 
one TLM mode and two TLH modes whose degenerate masses are, according to 
Eqs. (75]) and f90|). given as 

™? 2 = mL = 4 e 2 |0| 2 =4e 2 ||. (91) 

SU(3) case: 

Here 4> is diagonal in either one of the three SU(2) subgroups, and the gauge sym- 
metry is dynamically broken to U(l) 2 . We have 

A 3 

™ 2 -m? 9 =4e 2 



n2 — xll i2 



2vrT 



A 3 

m 2 3 = m 2 , = m 2 , = m 2 ., = e 2 



2ttT 



n3 — "H3 



A 3 

m T , = 4e 2 



2vrT 



2ttT 



'23 — xll 23 



A 3 

m 2 =4e 2 Ae 



m 2 3 = e 2 ^ • (92) 



2nT ' 

A| 
2ttT 

For T sufficiently far above T C: e we will see in Sec. 12. 3. 41 that e is practically inde- 
pendent of T. As a consequence, the TLH masses die off according to the power law 
in Eqs. (f9"Tj) and (f9"2"|) . Moreover, the energy density and the pressure of the ground 
state are only linear in T, P 9S = -p gs = -Au:K%T. We will show in Sec. 12.41 that 
the one-loop result for thermodynamical quantities is reliable on the 0.1% level 
throughout the electric phase. Thus the Stefan-Boltzmann limit P — j^T 4 " and 
p = ^T 4 (SU(2)) or P = ^T 4 and p = ^-T 4 (SU(3)) is reached very quickly 
apart from a factor arising from the extra polarizations of TLH modes. This result 
is in agreement with early lattice simulations using the differential method. 



2.3.2. Thermodynamical self consistency 

In this section we provide a conceptual basis for the notion of thermodynamical 
selfconsistency. 

As a result of the existence of a nontrivial macroscopic ground state, which is 
built of interacting calorons of topological charge-one, the ground-state physics and 
the properties of the excitations are temperature dependent. We have discussed in 
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Sec. l2.2l how the temperature dependence of the ground-state pressure and its energy 
density arises due to caloron interactions. These interactions are encoded in a pure- 
gauge configuration on the macroscopic level. For topologically trivial fluctuations 
Sa^ the two following properties are induced by the ground-state physics. First, in 
unitary gauge off-Cartan fluctuations are massive in a temperature dependent way, 
compare with Eqs. (|9Tj) and (|92]) . Notice that the quasiparticle masses are related 
to the ground-state pressure by their respective dependences on temperature if the 
temperature dependence of the effective gauge coupling e is known. Second, there 
are constraints for the maximal off-shcllncss of the fluctuations Sa^ arising from 
the compositeness of the ground-state field <p. Namely, a fluctuation Sa^, which was 
generated in a thermal equilibrium situation by the ground state, is not capable of 
destroying this ground state: In a physical gauge vacuum fluctuations or scatter- 
ing processes with momenta or momentum transfers larger than the (temperature 
dependent) scale \<f>\ are forbidden. 

Thermodynamical quantities such as the pressure, the energy density, or the en- 
tropy density are interrelated by Legendre transformations. These transformations 
can be derived from the partition function of the fundamental theory where the tem- 
perature dependence only occurs in an explicit way through the Boltzmann weight. 
A reformulation of the theory into a spatially coarse-grained Lagrangian, where cer- 
tain parameters are temperature dependent (implicit temperature dependences) by 
virtue of a separation into a ground state and (very weakly interacting) excitations, 
see Eq. (1731) , must not alter the Legendre transformations between thermodynamical 
quantities. For this to be true in the effective theory one needs to impose that in each 
transformation law the total derivatives with respect to temperature involves the 
explicit temperature dependences only. That is, derivatives with respect to implicit 
temperature dependences ought to cancel in a given Legendre transformation. 

A particular and essential Legendre transformation maps the total pressure onto 
the total energy density as 

P-rf-P. (93) 

If the effective theory has temperature dependent quasiparticle fluctuations of mass 
m a = c a m, where c a are dimensionless constants, and a ground-state pressure P gs , 
which can be regarded a fun ctio n of m, then the condition of thermodynamical 

I c 1 I 

selfconsistency is expressed as 

d m P = . (94) 

Eq. (|94|) assures that in Eq. ([93]) only derivatives with respect to the explicit tem- 
perature dependence of P contribute since jS; — drP + 9 m P4^. In d m P = 
the derivative of the pressure contribution arising from fluctuations cancels against 
that arising from the ground state. Eq. (|94|) governs the temperature evolution of 
the effective coupling e at any loop order that P is expanded in. 

The higher the loop order the more complicated the implementation of Eq. (|94D . 
In Sec. l2.3.4l we perform an analysis on the one- loop level (noninteracting quasipar- 
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Fig. 9. Diagrams contributing to the pressure when radiative corrections are ignored. A thick 
line corresponds to TLH modes and a thin one to TLM modes. The cross depicts the ground-state 
contribution that arises from interacting calorons. 

tides) which is more than sufficient for many practical purposes. A discussion of 
Eq. (|94[) on the two-loop level is given in Sec. 12. 4. 21 



2.3.3. Compositeness constraint and pressure at one loop 

We work in a physical gauge where TLM modes are transverse (two polarizations, 
Coulomb gauge with respect to the unbroken gauge group), and where TLH modes 
have three polarizations and do not interact with the pure-gauge configuration of 
the ground state (unitary gauge), for details see Sec. 12. 41 This is a physical gauge 
fixing which needs no introduction of additional fields since the Coulomb ghosts 
decouple from the dynamics. 

On the one-loop level, see Fig.[5J there are no interactions between the fluc- 
tuations 5a u . The only relevant compositeness constraint, related to the maximal 
off-shellness of a quantum fluctuation created by the ground state, thus is 

\p 2 — m 2 a \ < \<p\ 2 (Minkowskian signature) 
p 2 + ml < \cp\ 2 (Euclidean signature) . (95) 

SU(2) case: 

Before we discuss the thermal contribution to the one-loop pressure let us investigate 
what the constraints (|95p imply for the one-loop quantum correction — AVe- Setting 
m a = in (|95|) and considering two polarizations for TLM and three polarizations 
for TLH modes, an upper bound on |AVe| can be obtained as 



101 



\AV E \<\ ^ 3 l°g(rn)=^. ( 96 ) 



Thus we have 

AV E 



Ve 



101 V _ A s 



< 32tt 2 \Ae) ~ 32it 2 ■ (9?) 



Since > 13.867 when omitting AVe in the one-loop evolution of e, see Sec. 12. 3. 41 
this omission is justified: One-loop quantum corrections to the pressure are sup- 
pressed as compared to the tree- level result by a factor less than 2 x 10~ 6 . 
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Omitting the tiny quantum part, the one-loop expression for the pressure, in- 
cluding the ground-state contribution, reads: 

l4 



P(X E ) = -A%l—JL [2P( ) + 6P(2o)] + 2A £ 



(98) 



where A_e = ^f— and 

2T T V 2ttT 3 b v ; 

The (negative) function P(a) is defined as 

/•OO 

P(a)= / dxx 2 log[l - exp(-\/x 2 + a 2 )] • (100) 
Jo 

SU(3) case: 

Here a similar estimate for the quantum contribution to the one-loop pressure holds 
as in Eq. (|96p . and thus, again, we only have to consider the thermal part. It reads 

P(A E ) = - A e{^e t 4 ^(°) + 3 ( 4 ^(«) + 2 ^( 2 «))] + 2A ^} ( 101 ) 
with the same definitions as in the SU(2) case. 



2.3.4. One-loop evolution of effective gauge coupling 
Let us now implement the condition (|94|) . We have 

d m P = & 9 (oT) P = 0. (102) 

SU(2) case: 

Substituing Eq. (|98|) into Eq. (|102[) yields the following evolution equation 

_ 24A|a D(2a) 

Oa*E - ~ ,„ , 6 » ■ . 3 2 , (103) 

(27r) 6 1+ 24A^^ (2q) 



where the function D(a), see Fig-ITUl is defined as 

r°o 2 i 

D(a) = / <fa . (104) 

Jo vi 2 + a 2 exp(Va; 2 + a 2 ) — 1 

SU(3) case: 

The evolution equation is obtained by substituting Eq. (1101|) into Eq. Q102|) : 

^ = - » A V „° ( °' + 2D(2a) ■ (105) 

W l + i^(D(a) + 2D(2a)) 

Each of the Eqs. (|103|l and (|105| has two fixed points, one at a — and one at 
a = oo, see Fig.fTOl The points Xp.E = As(a — > 0) and A Ci £ = As(a = oo) are 
associated with the highest and the lowest temperatures, respectively, which are 
attainable in the electric phase. 
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Fig. 10. The functions aD(a) and a 2 D(a). 

Above \p,e the ground state would be trivial (topological fluctuations are ab- 
sent) and thus no tree-level quasiparticle mass were generated by a caloron induced 
Higgs mechanism. This re-introduces the problem of the diverging loop expansion 
as it is encountered in thermal perturbation theory and thus makes the ther- 
malized Yang-Mills theory inconsistent as an interacting field theory. We conclude 
that Xp^E marks the point in temperature where the field theoretic implementa- 
tion of the gauge principle breaks down. For a physics model, whose gauge group 
is a product of SU(2) and/or SU(3) groups, we expect that Xp t E ~ ^T^f~ wriere 
M P = 1.22 x 10 19 GeV is the Planck mass. 

As we shall see below, the point X c ,e, where each tree-level quasiparticle mass 
diverges, marks a transition to a (pre-confining) phase with condensed magnetic 
monopoles, dynamically broken dual gauge symmetries U{1)d (SU(2)) and U(l)|, 
(SU(3)), and isolated but instable center- vortex loops: The magnetic phase. 

Notice that the right-hand sides of Eqs. (|103|) and (| 1 05[) are negative definite. 
As a consequence, the solutions Ae(a) to Eqs. (|103p and (|105|) can be inverted, and, 
according to Eq. (|M)l . one obtains the evolution of the effective gauge coupling e 
with temperature as 

e(X E ) = ^a(X E )X 3 J 2 . (106) 

Z7T 

Inspecting the right-hand of Eqs. (|103p and (|105[) in the vicinity of the point X Cj e, 
it follows with Eq. (|106[) that e diverges logarithmically at X C: e- 

e(A B )~-log(A E -A C)B ). (107) 

In Fig.QT] a solution to Eq. (|103|) subject to the initial condition Xp iE = A^(a = 
0) = 10 is shown. We have noticed numerically that the low-temperature behavior 
of Xe(cl) is practically independent of the value Xp t E as long as Xp,E is sufficiently 
large. Let us show this analytically. For a sufficiently smaller than unity we may 
expand the right-hand side of Eq. (|103|) only taking the linear term in a into account. 
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Fig. 11. The solution Afi(a) to Eq. I|103|l subject to the initial condition Ap : g = 10 7 . 

In this regime, the inverse of the solution As (a) is of the following form 

"W 1 "(&)'■ (108) 

If Xe is sufficiently smaller than Xp e then this can be approximated as 

a oc X E 3/2 . (109) 

Thus the dependence in Eq. (|109|) is an attractor of the downward-in-temperature 
evolution as long as a remains sufficiently small: If we are only interested in the 
behavior of the theory not too far above A Cj b then it is irrelevant what the value 
of Xp t E is as long as it is sizably larger than A Cj £. This result is reminiscent of the 
ultraviolet-infrared perturbative decoupling property of the renormalizable, under- 
lying theory. Notice that the dependence of a on As in Eq. (| 109[) is canceled by the 
explcicit dependence of e on Xe in Eq. (| 1 06[) . Thus a plateau value e(Ae) = const 
is reached rapidly. 

In Fig. [12] the temperature dependence of e for SU(2) and SU(3), subject to 
the initial condition Xp t E — 10 7 , is shown for Xe < 500. Before we interpret our 
results a remark on the interpretation of the effective gauge coupling constant e is 
in order. Since e determines the strength of the interaction between topologically 
trivial gauge field fluctuations Sa^ and the macroscopic manifestation <f> of inter- 
acting calorons in the ground state e is not equal to the perturbatively generated 
gauge coupling constant g of the fundamental Yang-Mills theory. From Fig. [12] we 
see that the effective gauge coupling e evolves to values larger than unity. Naively, 
one would conclude that the theory is strongly coupled and that the one-loop evo- 
lution of e contradicts itself. This, however, is an incorrect conclusion. Because of 
compositeness constraints and the infrared regularization provided by TLH masses 
higher loop orders turn out to be very small as compared to the one-loop result for 
the pressure, see Sec. 12.4.21 
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Fig. 12. The temperature evolution of the gauge coupling e in the electric phase for SU(2) (grey 
line) and SU(3) (black line). The gauge coupling diverges logarithmically, e oc — log(As — \ c ,e), 
at X C:E = 13.867 (SU(2)) and \ CyE = 9.475 SU(3). The respective plateau values are e = 8.89 and 
e = 7.26. 



We interpret the fact that the gauge coupling constant e is constant for Xe siz- 
ably larger than X c .e as another indication for the existence of spatially isolated and 
conserved magnetic charges in the system. These charges are screened by calorons 
with a small holonomy, and thus e measures the effective magnetic charge of a BPS 
monopole which is given as g = — . The plateau values for e are e ~ 8.89 (SU(2)) 
and e ~ 7.26 (SU(3)). At X C , E = 13.867 (SU(2)) or X C , E = 9.475 (SU(3)) both the 
core size ~ - of a screened BPS monopole and its mass ~ — u max = ^7- vanish, 

e " e max ep ' 

see Sec. 13.1.11 (Notice that monopoles are very massive at high temperatures.) Thus 
monopoles are not well separated anymore at \ c ^e because it is extremely easy to 
move them. Local magnetic charge conservation is violated since in a typical vol- 
ume /3 3 the number of monopoles no longer is conserved. The increasing mobility of 
monopoles and the increasing violation of local charge conservation can be seen in 
the evolution of e for A# \ X c .e where an increase of e as compared to the plateau 
value is observed, see Fig. [12] and Eq. (|107|l . 

For Xe \ X Ct E TLH modes decouple (their masses diverge) and thus the (small) 
interaction between TLM modes dies off. This is the macroscopic manifestation of 
the fact that the magnetic charge of (dynamical and screened) BPS monopoles 
vanishes at A Cj £ making them unavailable as 'scattering centers' for TLM modes. 
This is, indeed, s een as a result of a two-loop calculation of the SU(2) pressure in 
the electric phase EBHH1 see Sec. l2A2l 
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Fig. 13. Diagrams for the TLM polarization tensor at one loop. Thick and thin lines denote TLH 
and TLM propagation, respectively. 



2.4. Radiative corrections for SU(2) 
2.4.1. Electric screening mass for TLM modes 

In this section we investigate for SU(2) the one-loop contribution to the electric 
screening mass for the TLM mode. 

In unitary gauge the analytically continued pro pagator of a free TLH mode 
Dj l ]f ab IJ,0 (k, T) is that of a massive vector boson E21 



j-vTLH.U.O/l rn\ r ( k^k^ \ 

D^ ab (k,T) = -S ab ^- im ,jx 

2^(fc 2 -4(eH)> s (|fc |/T) 



fc 2 — 4 (e|0|) 2 



, (110) 



where iib(i) = C J_ 1 is the Bose distribution and fco = ± y k 2 + 4 (e|</>|) 2 . The 
electric (or Debye) screening mass mn is related to the 00-component of the 
polarization tensor H )JjU (k) in the limit fco = 0, k — > 0. H fJ , ll (k) is calculated according 
to the diagrams in Fig.[13l The vertices for the interactions of TLH and TLM 
modes are the usual ones. In addition to the compositeness constraint (I95|) we have 
constraints associated with the maximally allowed momentum transfer in a fo ur 
vertex. Let the three independent momenta be p\,p2, and p 3 . Then we have^t 

I0l +P2?\ < |0| 2 , (s channel) \(p 3 - Pl ) 2 \ < |<?(>| 2 , (t channel) 

|(p 2 -p 3 ) 2 | < |0| 2 , (^channel). (Ill) 



One can easily see that these constraints together with the constraint (|95|) do not 
allow for any (neither from quantum nor from thermal propagation of the interme- 
diate states) contribution of the local (or tadpole) diagrams in Fig. [13] in the limit 
fco = 0, k — > and for e p i atcau = 8.89 (SU(2)). A calculation of the nonlocal diagram 
reveals that only thermal intermediate states contribute in the limit fco = and 
for e = 13.867. Notice that in perturbation theory the nonlocal diagram does not 
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Fig. 14. Nonvanishing two-loop diagrams contributing to the pressure. Thick lines denote prop- 
agators of TLH modes, thin lines those of TLM modes. 

contribute to U 00 (k = 0, k). The result for n o(fc = 0, k) read^l 



n 00 (fc = 0,fc) = 

2tt 1 2 \ \k\ + 4( e |0|) 2 



\k\ 



\kf 



(4(e|# 



2\2 




(112) 

Thus n o(fco — 0, k) is purely imaginary and diverges for \k\ — > 0: The electric 
screening mass mry, which is the positive square root of IIoo(fco = 0, fc — > 0), has 
an infinite and positive real part at finite coupling e: Static electric fields of long 
wavelength are completely screened by calorons. For e — > oo the Boltzmann factors 
in the integrals in Eq. (1112[) make IIoo(fco = 0, k) vanish at any momentum k. 



2.4.2. Two-loop result for the SU(2) pressure 

The nonvanishing two-loop diagrams contributing to the SU(2) pressure in the elec- 
trie phase llMJ are shown in Fig. [141 Because of the strong screening of near-to-static 
electric modes, compare with Sec. 12.4.11 the TLM propagator in Coulomb gauge can 
safely be approximated as 



D 



TLM.0/1 rri\ C 



P, 



I 



27TS(k 2 )n B (\k \/T) 



.UuU v 
1—*=; 

fc 2 



(113) 



e The author would like to thank Ulrich Herbst for performing this calculation. 
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where 

pT _ pT _ pT _ n 

= % ~ , (H4) 
k l 

k Q = ±\k\, and = S 0fl . 

The result of a calculation of the diagrams in Fig. [14] was published in SSI. We 
do only outline this (lengthy) calculation here. 

The following nomenclature is useful. Each diagram can be split into contribu- 
tions arising from the vacuum (v), the Coulomb (c) (oc u^u v in Eq. (|113[) ). and 
the thermal (t) parts of the involved propagators. Moreover, if a TLH propagator 
contributes to a given diagram then this situation is abreviated by H , in the other 
case by M. With e being larger than the one-loop plateau value e platoau = 8.89 the 
compositeness constraint in Eq. ([95]) allows for the five following two-loop radiative 
corrections to the pressure only: 

Up™, ^AP™, iAP™, \aP™, 1(AP™ + AP™). (115) 

The ratio of the two-loop corrections and the one-loop result (excluding the ground- 
state contribution) is plotted in Fig.[l5j For gp 1 — AP^ M the plot represents only 
an upper bound for the modulus of the correction, all other plots are exact results. 
The dominating correction 4p 1 - — AP™ M arises from the nonlocal diagram in 
Fig. 1141 It is negative. The dip is microscopically explained by the increasing effect 
of TLM modes scattering off of decreasingly massive magnetic monopoles close to 
the phase transition at X c ^e = 13.867. The effect of this scattering is suppressed 
for Xe S> X c .e since then the monopoles are dilute and massive scattering centers. 
(Recall that their mass is oc — after screening.) Notice, that the presence of massive 
but dilute scattering centers causes the contribution gp 1 - — AP^ M to remain finite 
but small for asymptotically high temperatures. 

A comment concerning thermodynamical selfconsistency on the two-loop level is 
in order. Recall that on the one-loop level we have obtained an evolution equation 
from the requirement of thermal selfconsistency d m P = 0. This gave a functional 
relation between temperature and mass which could be inverted for all tempera- 
tures in the electric phase. After the relation Eq. (|106|) between coupling constant 
e and mass a was exploited we obtained a functional dependence of the effective 
gauge coupling constant e on temperature. Equivalently, we could have demanded 
d e P = since e is the only variable parameter (apart from the scale Ae) of the 
effective theory in the electric phase. This would directly have generated an evolu- 
tion equation for temperature as a function of e. Because each diagram comes with 
a prefactor e 2 and due to the compositeness constraints radiative corrections AP 
to the pressure have a separate dependence on a and e, 



AP = T i AP(e,a,X E ), 



(116) 
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Fig. 15. Ratios (a) ^1—AP"" , (b) ^AP™, (c) > ^^AP™ , (d) > 
_l_ AP^ M , and (e) jp^- [ (AP^J IM + AP f f c SM ) as functions of temperature. 



where AP is a dimensionless function of its dimensionless arguments. To implement 
thermodynamical selfconsistency by demanding d m P — one has to express the 
explicitly appearing e in Eq. (|116p in terms of a by means of Eq. (|106|) and distinguish 
temperature dependences arising from a simple rescaling and those arising from the 
T dependent ground-state physics. For SU(2) we have 

p-'(«.W=T'^4»%. (117) 

The first factor on the right-hand sides of Eq. (|117[) arises from rescaling, so only 
the second factor needs to be differentiated: 

d a e 2 (a, \ E ) = 3ix (2a\ E + a 2 d a \ E ) • (118) 

After solving d m P — for the term d a X E we obtain a modified right-hand side of 
the evolution equation Eq. (|103[) . The two-loop evolution of e is work in progress^! 
The investigation of the screening masses and of the two-loop pressure for SU(3) is 
left for future work. The latter will be of a similar magnitude as in the SU(2) case. 
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3. The magnetic phase 

The electric phase is terminated at the temperature T c .e by the condensation of 
massless magnetic monopoles. These condensates can be described macroscopically 
by quantum mechanically and statistically inert complex scalar fields and pure 
gauges. The latter describe the interactions between monopoles which induce a 
negative pressure by the generation of isolated but collapsing magnetic flux loops 
(center-vortex loops). At T c _e one can not distinguish between the unbroken U(l) 
and the dual gauge group U(1)d (SU(2)) or the unbroken U(l) 2 and the dual 
gauge group U(l)|, (SU(3)): An exact electric-magnetic duality occurs. For SU(2) 
the point T Ct E is stabilized by a dip of the energy density. For SU(3) a similar 
stabilization takes place but with a much larger slope of the energy density: The 
phase transition appears to be weakly first order. It turns out that the electric Z2 
(SU(2)) or Z 3 (SU(3)) degeneracy of the ground state as it was observed in the 
electric phase is lifted: A unique ground state characterizes the magnetic phase. On 
the other hand, the expectation of the Polyakov loop, though strongly suppressed 
as compared to that in the electric phase, does not vanish entirely in the magnetic 
phase. This is a manifestation of the fact that the ground state in the magnetic phase 
does allow for the propagation of massive dual gauge modes despite the confinement 
of fundamental, fermionic, and heavy test charges by the monopole condensates. In 
that sense, the magnetic phase is only preconfining. 

The critical behavior in the vicinity of T c ^e is investigated and compared with 
results that seem to be related to the Yang-Mills theory by universality arguments. 

The monopole condensates are characterized by infinite correlation lengths. In a 
thermodynamical simulation on a finite-size lattice performed in the magnetic phase 
not much can be learned about thermodynamical quantities such as the pressure or 
the energy density which are sensitive to the infrared behavior of the theory. 

3.1. Prerequisites 
3.1.1. The BPS monopole 

Here we provide some facts about the 't Hooft-Polyakov monopole I 55 | 56 | - m |. ne gpg 
limit EQ] since we will need them in Sec. 13.1.21 We consider an SU(2) gauge model 
with the Lagrangian of Eq. ([75)) with the modification that the potential is absent. 
The BPS monopole is a static, particle-like solution to the equations of motion of 
this model saturating the BPS bound on the mass M. When centered at x = it 
is given as 

<^ = z a |0|F(p), oS = 0, a? = G(p) , (119) 

where r = \x\, Xi = —. \(f>\ = V '<j) a <j) a (r — > 00), and p = e\<p\r. The antimonopole 
solution is obtained by letting x — > —i in Eq. (|119[) . The functions F and G can be 
determined analytically. They are given as ^01 

F(p) = coth p-~, G{p) = 1 - . (120) 

p sinn p 
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The mass M of a BPS monopole or antimonopole calculates as 

M=—\<f>\. (121) 

e 

(Notice that in Eq. (|12ip \cf>\ is replaced by u ~ u max = 5 for a screened BPS 
monopole generated by the dissociation of a large-holonomy caloron, compare with 
Eq. I|13p.) The magnetic charge g is obtained by integrating the 't Hooft tensor 

= - dvfraft - -e abc rd^ b 8j c , {#> = , (122) 

e 

over a two-sphere 5*2 with infinite radius which is centered at x = 0. It is given as 

47T 

dE^F^ = ±— . (123) 

S 2 M=oo e 

In Eq. (|123[) denotes the differential surface element. 

In unitary gauge, where the color orientation of <j> a is 'combed' into a fixed 
direction in adjoint color space, the magnetic field Bi = T~§k associated with 
Eq. (|122p is accompanied by a Dirac string along this direction in position space: 
If the monopole lies inside an 5 2 then the magnetic flux through this S 2 of the 
Dirac string precisely cancels that of the hedge-hog magnetic field. If a monopole or 
antimonopole lies outside of an S 2 with infinite radius an finite distance b away from 
its surface and the monopole's or the antimonopole's Dirac string does not pierce 
the surface then the magnetic flux F through the surface, see Fig. [HI calculates as 

F = J dYi^ v T^v 

Y />2tt />oo 



2 />Z7T />00 

±— / d/3 dx \x\ cos a(x, b) 

e Jo J -00 



x 2 + b 2 



e J {x 2 + b 2 fl 2 

4-7T 

= ±— . (124) 
e 

We conclude that static monopoles or antimonopoles lying inside an S2 of infinite 
radius do not contribute to the flux through S% while they may contribute to the 
flux when situated outside of this S 2 - 



3.1.2. Derivation of the phases of macroscopic complex scalar fields 

There is one species of magnetic monopoles in the SU(2) case while there are two 
independent species of magnetic monopoles for SU(3). 

SU(2) case: 

In unitary gauge we consider an isolated system of a monopole and an antimonopole, 
which both are at rest and do not interact, outside of an S2 with infinite radius. 
We characterize their Dirac strings by unit vectors x m and £ a which point away 
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outside 



b 




P x 



inside 



Fig. 16. A BPS monopole in unitary gauge outside of an S2 with infinite radius. The Dirac string 
does not pierce the surface of the S2- 

from the core of the monopole and the antimonopole, respectively. Let P be the 
plane perpendicular to S2 such that the intersection line L = P n S2 coincides with 
the intersection line of S2 with the plane spanned by x m and x a . (The case where 
x m and x a lie in S2 is inessential for what follows.) Whether or not this system 
contributes to the magnetic flux through S2 depends on the angle <5 = /.(x m ,x a ) 
and on the angle 7 which the projection of x m onto P forms with L. 

A magnetic flux through S2 is generated if and only if either x m or x a alone 
pierces S2. For a given angle 6 the angle 7 is uniformly distributed. In the absence 
of a heat bath the probability of measuring a flux yora flux — ^ through S2 thus 
is given as j- . We conclude that for a given angle S the average plus or minus flux 
through S2 reads 



Notice that 5 and 2ir — S generate the same average flux F± , thus the restriction 
< 5 < 7T in Eq. (fl25|) . 

So far we have discussed the flux through S2 which is generated by an isolated 
monopole-antimonopole system with no interactions. To derive the phase of the 
macroscopic complex field tp describing the Bose condensate of such systems we cou- 
ple the system to the heat bath, project onto zero-momentum states (condensate) 
of the monopole-antimonopole system such that each constituent does not carry 
momentum and perform the massless limit e —* 00 which takes place for T < T c ^ 1 
compare with Eq. (|107p . (The rare case of a zero-momentum state with opposite 
and finite momenta of the constituents generates a closed and instable magnetic 
flux line, see Fig.[TJ This situation takes place if a large number of large-holonomy 
calorons dissociate into monopole-antimonopole pairs almost simultaneously inside 
a small spatial volume. On the macroscopic level, the thermal average over these 
flux loops, which collapse as soon as they are created, will later be described by a 
pure gauge configuration.) 



F± 




(0 < 5 < tt) . 



(125) 
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The monopole and antimonopole are generated by the dissociation of a large- 
holonomy caloron. According to Eq. (|25[) the sum of monopole and antimonopole 
mass, M m -\- a , is, after screening, given as 



M„ 



8tt2 

ej3 



(126) 



The thermally averaged flux of the zero-momentum system at finite coupling e is 
obtained as 



F±, th (6) = 4n J d 3 p6 (3 \p)n B (p\E(p)\)F ± 



± — / 



exp 



. p2 



(127) 



After setting p = in [exp Py M% +b + p 2 — lj and by appealing to Eq. (|126[) . 
the expansion of this term reads 



lim exp 



8tt* 
e 



18tt2 
2~e~ 



(128) 



Appealing to Eq. (|128[) . the limit e 
We have 



oo can now safely be performed in Eq. (|127|l . 



lim F± tth {5)=±- 



(0 < <5 < vr) . 



(129) 



The right-hand side of Eq. p29p defines the argument of the complex and periodic 
function / with 



/(- 



\tp\ IT 



(130) 



where C is an undetermined (complex) constant. (Recall, that 5 is an angle). Since 
/'s argument was obtained by a projection onto zero spatial momentum the only 
admissible nontrivial periodic dependence is that on the Euclidean time r. Without 
restriction of generality we can thus set £ = 5. Since / is periodic, /(r = 0) = 
/(t = 0), it can be expanded into a Fourier series: 



T 

27r«n— 



(131) 



where /„ are (complex) constants. According to Eq. (|130p // = |C| 2 is constant in 
r. Thus the only possibility in Eq. (|131[) is /„ = C <5 rml or /„ = C <5(_ m ) n for a fixed 
value of m. Moreover, only m = 1 is allowed since the physical situation generating 
the continuous parameter 8 does not repeat itself for 0<<5<^,^<<5<^,--- 
if m > 1 because no higher-charge monopoles exist. (Recall that only calorons with 
a large-holonomy and topological charge unity are allowed to contribute to the 
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thermodynamics in the electric phase.) We conclude that the equation of motion 
satisfied by / is: 

3?/ <M¥) 2/( i ,=0 - (132) 

3.1.3. Derivation of the modulus of macroscopic complex scalar fields 

What about c^'s modulus? The reasoning is completely analoguous to that for the 
derivation of 0's modulus. First, since ip is composed of massless, noninteracting 
monopoles being at rest its r dependence ought to be BPS saturated. Second, for 
f to have the phase exp[±27ri^] the right-hand side of its BPS equation ought to 
be linear in ip. Third, for the existence of smooth deformations (3 — > (3 + 8(3, subject 
to a perturbative expansion in ^ away from a phase boundary, the right-hand side 
of the BPS equation must be analytic. Fourth, we assume that a scale Am is given 
externally. Fifth, no explicit dependence on (3 may appear on the right-hand side 
of the BPS equation since at zero momentum the temperature dependence of the 
mass of a monopole cancels against the factor (3 in the Boltzmann weight. The only 
possibility for the BPS equation satisfying these conditions is 

o T <p = ±i i i, = ±1—^ . (133) 

wr v 

From Eq. (|133|) it follows that ip's modulus is given as 



M = V^- (134) 

The right-hand side of Eq. (|133|) defines the square root of the potential Vm- In 
the absence of interactions between (screened) monopoles the effective theory for ip 
thus reads 

S v = [ dr [ o?x(\l\^d T ip+\^] . (135) 



o 



2 r r 2 if ip 



SU(3) case: 

Here we have two independent monopole species which do not interact. The situation 
for each species is completely analogous to the SU(2) case. The two macroscopic 
fields if i and tp% both satisfy the BPS equation (|133|) and are given as 



<Pi(r) = <P2(t) = \/ exp 



(136) 



In the absence of interactions between monopoles the effective theory for <fi,f>2 
reads 

r /i i 1 A 6 1 A 6 \ 



'0 



S v = dr d A x[- d^dr^ + - d T ip 2 d T ip 2 + - + - . (137) 

' ' \2 2 2 ip x ipx 2 ip 2 f>2j 
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3.2. A macroscopic ground state 

In close analogy to the electric phase we now derive the full macroscopic ground- 
state dynamics. First, we establish the quantum mechanical and statistical inertness 
of the fields ip (SU(2)) and <^i,<£>2 (SU(3)). Subsequently, we solve the equations of 
motion for the topologically trivial dual gauge field in these backgrounds to obtain 
pure-gauge configurations describing, on a macroscopic level, the interaction be- 
tween monopoles in the ground state. It will turn out that the Polyakov loops when 
evaluated on these pure-gauge configurations are trivial: The electric Z 2 (SU(2)) 
and Z3 (SU(3)) degeneracies observed in the electric phase no longer exist. Thus 
the magnetic phase confines fundamental test charges. 

The ratios of the mass-squared of potential ip- field fluctuations with \tp\ 2 and T 2 

are 

where X M = ^Z. We will show below that X M > 7.075 (SU(2)) and A M > 6.467 
(SU(3)). Thus ip-&e\d (SU(2)) and ipi,tp 2 -&eld (SU(3)) fluctuations neither exist 
on-shell nor off-shell. 



3.2.1. Pure-gauge configurations 
SU(2) case: 

The topologically trivial sector is coupled to ip in a minimal fashion, and the fol- 
lowing effective action arises 

S = J dr J d 3 x 

where = d u a® — d v a® denotes the Abelian field strength of the dual gauge field 
a® and = + ig a® denotes the covariant derivative involving the magnetic 
gauge coupling g. 

Since the field ip does not fluctuate it is a background to the macroscopic gauge- 
field equations of motion which follows from Eq. (I139P : 

= *<? [V^P<P " &>M • (140) 
A pure-gauge solution to Eq. (|140p with D u tp = is given as 

a^=±S^. (141) 

On ip and on a^' bg only the potential does not vanish in Eq. (|139[) : Interactions 
between magnetic monopoles create a nonvanishing energy density p 9S and pressure 
P gs where 

p 9S = nA 3 M T = -P 3S . (142) 



(139) 
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We shall see in Sec. 14.1.11 compare with Eq. (|165|) . that the negative ground-state 
pressure in Eq. (|142[) originates from center-vortex loops which collapse as soon as 
they are created. 
SU(3) case: 

The situation is the same as for SU(2) except that we have gauge dynamics subject 
to U(l)|) and not only U(1)d. The effective action reads 



S = f dT I d 3 " 







1 A?, 1 A?, I 

(143) 



!^ ~ !^ ^ 1 A|, f 1 A?, 



2 p ' iri p ' iri 2 2<p m 2^2 

The Abelian field strengths l5 G^„ 2 and the covariant derivatives T> fij i,T> flj 2 are 
defined as for the SU(2) case with the replacements a® — ► ar' 1) gc — > a^ 2 , respec- 
tively. (The magnetic gauge coupling g is universal since both species of monopoles 
couple with the same strength to their respective gauge field.) 

The equations of motion for the fields o>ui> a u,2 m * ne background of the fields 
tpi, <p 2 are 

d ^G^ v l = ig [D v ,iipi<p - ^l^i^i] , ^G^ 2 = ig [V v ,2<p2<P ~ ^2^,2^2] ■ (144) 
Pure-gauge solutions to these equations with D v iipi = D V 2 i P2 = are given as 

On ipi, tf2 and on a®f 9 , a®'^ 9 only the potentials do not vanish in Eq. () 143[) : Inter- 
actions between magnetic monopoles create a nonvanishing energy density p 9S and 
pressure P 9S where 

p 9S =2ttA 3 m T= -P 9S . (146) 

Again, the negative ground-state pressure in Eq. (|146[) originates from center- vortex 
loops which collapse as soon as they are created. 



3.2.2. Polyakov loop and rotation to unitary gauge 
SU(2) case: 

ig /,f dr af 



The Polyakov loop P = exp 
figuration in Eq. (|141[) . reads 



, when evaluated on the pure-gauge con- 



P = cxp[±ig 



\ 27r l 1 

gp 



(147) 



A gauge rotation a^ bg 



(dfiSV) to unitary gauge ip 



is mediated by the U(l) group element Q = exp 



The Polyakov loop 



P is invariant under this gauge transformation. We conclude that the electric Z2 
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ground-state degeneracy, which was observed in the electric phase, no longer exists 
in the magnetic phase: The ground state confines fundamentally charged, heavy and 
fermionic test charges. 

SU(3) case: 

Here the Polyakov loop P is a product of the Polyakov loops Pi and P2 computed 
on the respective pure-gauge configurations a®'^ 9 and a^'^ 9 in Eq. (|145[) . We have 



'Jo dT ^ ] 



P = P1P2 = exp[±2i 3 / dr — } = 1 or P = PiP 2 = exp[0] = 1 . (148) 
Gauge rotations to unitary gauge ip\ — \ip\\ — ip2,a^\ g = a 1 ?^ 3 — are mediated 



by the U(l) group elements fli = cxp 



±27U-§ 



= Q?. Again, P is invariant under 



these gauge rotations. We conclude that the electric Z% ground-state degeneracy 
does not exist in the magnetic phase: Fundamentally charged, heavy and fermionic 
test charges are confined by the monopole condensates. 



3.3. Excitations 

3.3.1. Mass spectrum of thermal quasiparticles 

The dual Abelian Higgs mechanism generates tree-level quasiparticle masses m and 
mi,m 2 for the fluctuations Sa^ (SU(2)) and Sa^^Sa^ (SU(3)), respectively. We 
have 

m = g\ip\ = TOj = 77i2 = clT , (a = 2ir g \^f^ 2 ) (149) 
where A M = 1 T77- 



3.3.2. Thermodynamical self consistency and evolution equation 



Due to the absence of interactions between the dual gauge fields the thermodynamics 
of the magnetic phase is exact on the one-loop level. Again, a magnetic modification 
of the compositeness condition Eq. (|95|) applies. 

Let us first compare the contribution AVm of quantum fluctuations to the pres- 
sure arising from dual gauge modes with the tree- level result —1/2 Vm 
(SU(2)) and —1/2 Vm = ~2n A 3 M T (SU(3)). In both cases we have 



nA M T 



M 



X" 3 



V M - 24^ ■ ( 15 °) 

Considering that A M > 7.075 (SU(2)) and A M > 6.467 (SU(3)) this is smaller 
than 1.2 x 10 -5 and 1.6 x 10 -5 , respectively. Thus the quantum contribution to the 
one-loop pressure can safely be neglected. 

For SU(2) the thermal contribution to the pressure reads 

(2tt)6 ■ W+ 2 _ 



P(A M ) = -A 



(151) 



1, 2008 20:36 WSPC/INSTRUCTION FILE F 



56 Half Hofmann 




Fig. 17. The evolution of the effective gauge coupling g in the magnetic phase for SU(2) (thick 
grey line) and SU(3) (thick black line). At A CjM = 7.075 (SU(2)) and X C>M = 6.467 (SU(3)) g 
diverges logarithmically, g ~ — log(Ajif — \ c ,m)- 



where the (negative) function P(a) is defined in Eq. (|100|) . The SU(3) pressure is 
just twice the SU(2) pressure. From the condition d a P — of thermodynamical 
selfconsistency the following evolution equation arises for both SU(2) and SU(3): 



Oa^M — —- 



D(a) 



(27T) 



1 



(2tt)6 



■25(a) 



(152) 



where the (positive) function D(a) is defined in Eq. (|104[) . In analogy to the electric 
phase the Am dependence of the gauge coupling constants g is obtained by inverting 
the solutions to Eq. (|152[) and by subsequently using the relation between g, Am and 
a in Eq. (] 149[) . The temperature evolution of g is shown in Fig. 1171 



3.3.3. Interpretation of results 

The magnetic gauge coupling g increases continuously from g = at the electric- 
magnetic phase boundary (T = T c< e) until it diverges logarithmically at T Cj m- A 
variation of the magnetic coupling with temperature is not in contradiction with 
magnetic charge conservation since no isolated magnetic charges appear in the mag- 
netic phase: Magnetic monopoles either are condensed or they conspire to form in- 
stable magnetic flux loops (center- vortex loops), see Fig.[TJ From Fig. [IT] one can see 
that the magnetic phase is more narrow for SU(3) than it is for SU(2). Despite the 
fact that the ground-state degeneracies with respect to the electric Zi symmetry 
(SU(2)) and the electric Z% symmetry (SU(3)) are absent in the magnetic phase the 
fully averaged Polyakov loop (including the massive dual gauge-mode excitations) 
does not vanish at T c ,e, see Sec. 13. 41 The expectation of the Polyakov loops only 
vanishes at T Ci m where all dual gauge modes are decoupled because of a diverging 
mass. If we take the mass of the dual gauge modes to be the order parameter for the 
dynamical breaking of U(l)u (SU(2)) and U(l)fj (SU(3)) then the electric-magnetic 
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transition is second order with mean-field exponents in both cases. The best one can 
do to relate this order-parameter to electric or Z3 restoration is to look for the 
point where its exponent is least sensitive to the length AAm of the fitting interval. 
The expectation of the 't Hooft loop, which is a dual order parameter for complete 
confinement and which is nonzero if center-vortices are condensed, vanishes inside 
the magnetic phase where center- vortex loops are isolated and instable defects. At 
T c m the expectation of the 't Hooft loop jumps to a finite value. The associated 
transition is, however, neither second nor first order but of the Hagedorn type see 
Sec. 14.41 Because of the infinite correlation length (massless condensed magnetic 
monopoles) a finite-size lattice simulation of the order parameter as well as infrared 
sensitive thermodynamical quantities such as the pressure is problematic. 



3.4. Polyakov loop in the electric and the magnetic phase 

In this section we show, on a qualitative level, that the expectation of the Polyakov 
loop (P), which is an order parameter for the confinement-deconfinement transition, 
is finite both in the electric phase and the magnetic phase. Deep in the magnetic 
phase (P) is, however, strongly suppressed as compared to its value at T Ct E- 

In each phase the Polyakov-loop expectation (P) can be computed in unitary(- 
Coulomb) gauge. Let us first discuss the electric phase. The sector Sj^e in the 
effective action (|T3|) - which involves fluctuating fields, is given as 



S f,E - 



drcPa 



(153) 



Here a = 

mi — 777-2 



1, 2, 3 and mi 



) m 3 



T m 4 



TO2 

\m 5 



> 0,7773 = for SU(2) and a = 1, • ■ ■ ,8 and 
= \m§ > 0,7777 = ma = for SU(3). Since in 



unitary gauge the Polyakov loop is unity in the ground state no direct ground-state 



contribution arises in the expectation (P): 
cancels that in the denominator. We have 



— Z f,E x 



The associated factor in the numerator 







j [ddafj] exp 


ie / dr Sa4 




Jo 



exp[-5/,B] 



(154) 



where [d5a^] denotes the path-integral measure and Zf t E = / [dSa^] exp[— 

Since the fluctuations <5a M are periodic in r they can be decomposed into a 
Matsubara sum: 



6a M (r,f) 



n— 00 

E 

n— — cx 



exp 



27T777- 



0. 



5a^ n {x) 



(155) 



Modes with 77 ^ in Eq. (|155[) render the Polyakov-loop phase in Eq. (|154[) to be 
unity and are action-suppressed. Zero modes (77 = 0) contribute to (P) sizably if 
they are not action-suppressed. This is the case if and only if both of the following 
conditions are met: (i) as compared to T some or all masses m a in Eq. (j 1 53(1 are 
small and (ii) di5a^x>{x) is small compared with T 2 and the field configuration is 
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still localized in space. Here d% denotes a spatial derivative. In the electric phase 
there are massless modes, the conditions (i) and (ii) can be satisfied, and thus a 
finite Polyakov-loop expectation emerges. 

A similar consideration can be performed for the magnetic phase. Deep inside 
this phase condition (i) is badly violated since the mass of all dual gauge modes is 
much larger than T by virtue of Fig.fTTl and the mass formula in Eq. (] 149|) . Thus 
(P), though nonvanishing, is strongly suppressed deep inside the magnetic phase 
as compared to its value at T c ^e- For quantitative results the average (P) can be 
performed on a lattice or analytically by using the respective effective theory for 
the electric and the magnetic phase. 

3.5. Critical behavior at the electric-magnetic transition 

The electric-magnetic transition, which goes with the dynamical breakdown of 
U(l)r> (SU(2)) and U(l)|, (SU(3)), is second order in both cases. The difference 
is that in the SU(3) case the magnetic phase is more narrow than for SU(2), see 
Fig.1171 and that the peak-value of the specific heat per volume is much larger for 
SU(3) than it is for SU(2), see Fig. [29] In addition, the entropy density, which mea- 
sures the mobility of dual gauge modes and which is used to extract an apparent 
latent heat on the lattice drops much more rapidly in the magnetic phase for 
SU(3) as compared to the SU(2) case. This is a plausible explanation for the ap- 
parent first-order nature of the confinement-deconfinement transition observed in 

loo QO 

lattice simulations ° [ 

The order parameter for the electric-magnetic transition is the mass oc osAm 
of the dual gauge bosons. (The monopole mass vanishes like an inverse logarithm 
on the electric side of the transition, and thus it is not an order parameter.) The 
following model applies to the behavior of a Am close to the critical temperature 
A c = 9.24 (SU(2)) and A c = 6.81 (SU(3)): 

a\ M (\M) = K\\ M -Xc\ v , (156) 

where K and v are constants, and A c is the critical temperature T c< e in units of -^ d -. 
By demanding continuity of the pressure across the electric-magnetic transition, see 
Sec.El we derive A c = 8.478 (SU(2)) and A c = 7.376 (SU(3)). 

The magnetic phase is not entirely confining (dual gauge bosons propagate al- 
though fundamental, heavy, and fermionic test charges are confined), and thus the 
expectation of the Polyakov loop (P) is not exactly zero. The best one can do in 
order to compare t he beha vior of a Am to the behavior of (P) inferred from uni- 
versality arguments I 69 | 70 | j g ^ Q Iqq^ f or the point where the fitted value of v in 
Eq. (|156p is least sensitive to a variation of the length AAm of the fitting interval. 

To perform the fit to the model in Eq. (|156p we have used Mathematica's Non- 
linearFit function which is contained in the statistics package. The function o(Am), 
subject to the initial conditions a(\ c = 8.478) = (SU(2)) and a(\ c = 7.376) = 
(SU(3)), was generated by an inversion of the corresponding numerical solutions to 
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AA M 

Fig. 18. The critical exponent u for the mass of the dual gauge modes and for the magnetic- 
electric transition as a function of the length A Am of the fitting interval. The left panel is the 
SU(2) and the right panel the SU(3) result. 

Eq. p52[) . (A step-size 8a = 5 x 10~ 9 was used in the Runge-Kutta algorithrrQ. ) In 
Fig. HH the fitted exponent v is shown as a function of AAm- Two things are im- 
portant to observe. First, the magnetic-electric transition is second order with the 
mean- field exponent v — 0.5 for both SU(2) and SU(3). Second, for each case there 
is a point AA^ of least sensitivity for v under variations in AAa/. For SU(2) we have 
AA* f = 0.28 ±0.03 and z^(AA^ f ) = 0.576 ±0.008, and for SU(3) AA M = 0.16 ±0.02 
and K AA iW) = 0.572 ± 0.006. By universality we expect f(AA M ) for SU(2) to be 
close t o the exponent u IM for the corresponding order parameter of a 3D Ising model 
l 69 l 7 °l . One has v lM ~ 0.63. The SU(2) exponent i/(AA M ) only deviates by about 
8.5% from v lm . 

4. The center phase 

4.1. Prerequisites 

4.1.1. ANO vortex in the BPS limit 

Just like the isolated defects in the electric phase are screened BPS monopoles, 
the isolated defects in the magnetic phase are screened and closed magnetic flux 
lines (vortex- loops) . In Fig.[TJ see also ^Bl^ we have given a figurative interpretation 
of these flux lines: They are composed of magnetic monopoles and antimonopoles 
which move into opposite directions. Thus there is a net magnetic current in the 
vortex core. A vortex core can be viewed as locations in space where U(1)_d (SU(2)) 
or one of the factors in U(l)|, (SU(3)) are restored. Hence the picture of isolated 
monopoles, contributing to the magnetic current, applies to the vortex core. 

We have explained in Sec.[T] why the magnetic flux carried by a vortex-loop is 
independent of the state of motion of a particular monopole contributing to the vor- 
tex: The amount of flux carried by a vortex solely is a function of the charge of each 
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BPS monopole contributing to it. Because large-holonomy calorons of topological 
charge one and thus monopoles with one unit of magnetic charge only are thermo- 
dynamically excited Abrikosov-Nielsen-Olesen (ANO) vortices are center vortices 
in the magnetic phase of an SU(2) or an SU(3) Yang-Mills theory. 

A mesoscopic description of a static ANO vortex in the BPS limit is given by 
the action Eq. (| 1 39[) when omitting the potential for ip. (This potential measures the 
energy density of the macroscopic ground state and thus must be subtracted when 
discussing the typical energy of a solitonic configuration on a mesoscopic level.) 
The following cylindrically symmetric (with axis along the x% d irection) and static 
ansatz for the gauge field a D is made to describe the vortex 

of = , aP = Cijkfj e k A{r) , (157) 

where f is a radial unit vector in the X1X2 plane, and e is a unit vector along the 
X3 direction. Writing tp = |y|(r) exp[z#], the equations of motion for |y|(r) and A(r) 
read 

-HH^+Q-'tf™- ' (i58) 

-i(-T^) + ^\ 3 (9^--)=Q- ( 159 ) 

dr \r dr I \ r I 



We keep in mind that \<p\(r — ► 00) = \J see Eq. (|134| ). Let us first search 

for BPS saturated solutions to the system (|158I159[) . The question is under what 
condition a solution to the first-order equation 

l M =(i- 5 A)M (160) 

also solves Eq. (1158[) . Substituting Eq. (|160[) into Eq. (|158|) , we observe that 

A=-r4-A. (161) 

dr 

The solution to Eq. p^T|) is 

, , , const , , 

A(r) = . (162) 

r 

Substituting Eq. (| 162|) into Eq. (| 1 59[) , the constant in Eq. (| 162[) is determined to be 
-. Eq. (fT60|) is solved for r > by 



W\{r) = s]^-. (163) 

We have found an analytical solution to the system (|158I159[) for r > which has one 
unit of magnetic flux F v (r) = ^ = § c dz^ a D , where C is a circular curve of radius 
r in the X1X2 plane centered at x\ = x-i = 0, and which has a vanishing vortex core. 
The energy density p v {r), when evaluated on the configuration A(r) = j^, |y?|(r) = 

/ A3 ~Q 1 , 

y 2k ; reduces to that of the magnetic field H(r) — ^( r ) : (^y Stoke's 
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theorem the magnetic field H must be proportional to a two-dimensional (5-function 
at r = 0. Thus the energy per unit vortex length diverges on the configuration (|162[) 
and ([163]).) 

Pv (r) = ^H 2 (r) = 0, (r>0). (164) 
The (isotropic in the X1X2 plane) pressure Pv{f) outside the vortex core is given as 



p »W = -^ti. ( r >°)- ( l65 ) 

2 2tt g A r z 



Eq. (|165p is the mesoscopic reason for the macroscopic results in Eqs. (|142|) and 
(|146p . Because of the negative pressure in Eq. (|165|) vortex-ioops start to collapse as 
soon as they are created at finite coupling g. (The pressure is more negative inside 
than outside of the vortex-loop.) Notice that in the limit g — > 00 we have P v (r) — > 0: 
For temperatures below T c .m vortex- loops do exist as particle- like excitations. 

4.1.2. Leaving the BPS limit 

Let us now discuss how the solutions in Eqs. (|163jl and (|162p are deformed when 
the BPS limit is left at finite coupling g. In this case only approximate analytical 
solutions to the second-order system (|158[) and (|159[) are known 1^3. Assuming \tp\ 
to be constant, which is viable sufficiently far away from the core around r = 0, the 
solution to Eq. (|159[) reads 



A{r) = — - \<p\ K x {g\ip\r) -> — - [yU/ f" , exp[-g|y]r] , (r ->• oo) , (166) 

where ifi is a modified Bessel function. (Notice that the 1/r divergence at r = 
of the solution in Eq. (|162[) is absent in the configuration in Eq. (|166[) .) Now \<p\ is 
not constant inside the vortex core but smoothly approaches zero for r — > 0. So 
there is a gradient contribution from \ip\ to the energy per unit length along 
the vortex where R denote the typical radius of a vortex loop. Let us first calculate 
the magnetic energy per unit length ■ One has^^ 

l^=2j dr ^rH 2 (r)=Tr\v\ 2 ^ dyKl{y)y = -\ V \\ (167) 

where \tp\ is given in Eq. (|163|) . The gradient contribution -§^jjr is comparable to 

E, 



2nR 

2^W' Thus the typical energy E v of the vortex loop is obtained by multiplying 



E m v , e. 



+ with the typical circumference 2nR ~ -r-r of the loop. We have 



2-kR ' 2ttR " v-ii^uiiu^^i, 

^~2f M 2 X^=,M. ( i 68) 

From Eq. I|168p we conclude that vortex loops become massless in the limit g — > 
oo. For r S> the (isotropic in the X1X2 plane) pressure P v (r) of the vortex 
configuration is still given by Eq. (|165[) : For finite coupling g vortex loops collapse 
as soon as they are created. 
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Fig. 19. Two center-vortex loops of opposite flux being pierced by an S± of infinite radius. 



4.2. Derivation of the phase of a macroscopic complex scalar field 

We consider a pair of center- vortex loops that are pierced by a circular contour C of 
infinite radius, see Fig.[l9j The total center flux F± o through the minimal surface 
spanned by C is 

f ±2ZL 

*±,o - < Q 9 (169) 

depending on whether at finite coupling g the loop A (B) collapses to nothing well 
before the loop B (A) or whether this roughly happens at the same time. (In any 
case, vortex loops which start out without getting pierced by C do not contribute 
a center flux through C.) 

Unlike in the case of a pair of a monopole and an antimonopole discussed in 
Sec. l3.1.2l the flux F±.o takes discrete values. In analogy to the derivation of a macro- 
scopic monopole condensate we may investigate the thermally averaged flux of the 
vortex-antivortex (spin-0) system in the limit where there is no spatial momentum 
of this system and where g — > oo: 

km F ±i0 . tb (8) = 4tt / d 3 p5W{p) n B (f3\2E v (£)\) F ±fi 

9 — oo J 

3 /2 

= 0,±^ T = I ±^£. (170) 
The phase of a macroscopic complex field $ is defined as 

r^Ca:) = lim ( exp[ig <f dz M (a D r]\ (171) 

W 9^°° \ Jc(x) I 

where T is a complex constant, and (a D )^ denotes the gauge field of a center vortex. 
The expectation on the right-hand side of Eq. (|171[) is proportional to the expecta- 
tion of 't Hooft's loop operator (Green functions of this operator change their 
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Z 2 monopole 




Fig. 20. The creation of an isolated Z2 monopole by self-intersection of a center-vortex loop. 

phase by —1 (SU(2)) and exp[±^p] (SU(3)) under an exchange of the order of 
any two of their arguments: A feature which is familiar from Green functions of 
fermionic fields.) The possible values of $'s phase are parametrized by the average 
center flux limg^oo F± o- th (5) in Ea.([l70ll. 

According to Eq. (|170p the condensate $ of center- vortex loops is determined by 
discrete parameter values which can be normalized as f = ±1, 0. Recall that at A Cj m 
vortex loops start to be stable excitations since their pressure P v (r) is zero outside 
the (infinitely thin) vortex core. Once the field <E> acquires a nonzero modulus its 
phase is observed to jump locally in space. (Each jump corresponds to a stable 
vortex loop travelling in from infinity and getting pierced by C) Thus we are led to 
interpret jumps in $'s phase as creation processes for (fermionic) particles. (There 
are two degenerate polarizations of these particles: The two possible directions of 
center flux in a given vortex-loop. By travelling in from infinity the vortex loop 
makes $'s phase jump twice: A created unit of flux is associated with a forward 
jump (piercing by C) while a backward jump corresponds to minus this flux (no 
piercing by C, center- vortex loop lies inside C).) If a single center- vortex loop is 
created with sufficiently large momentum then a part of its kinetic energy can 
be converted into the mass of self-intersection points by subsequent twisting. Self- 
intersection points are Z2 magnetic monopoles, each contributing ~ Ac to the mass 
of the state where Ac is a mass scale. Twisting does not alter the fact that only two 
possible polarizations (spin-1/2 fermions) occur. (The magnetic flux is reversed by a 
Z2 monopole l-^H see Fig.[20l) We conclude that the mass spectrum m n of fermionic 
excitations is equidistant: 

m„~nAc, (n = 0,l,2,---). (172) 

The process of fermion creation violates the spatial homogeneity of the system 
and thus thermal equilibrium. Fermion creation, that is, the process of sucking in 
stable center- vortex loops from infinity, can only go on so long as the energy density 
provided by the field $ is finite. Thus the field $ must eventually relax to one of the 
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possible zero-energy minima of its potential. This phenomenon is generally known 
as tachyonic pre- and re-heating 1241 
SU(2) case: 

The symmetry, which is dynamically broken by center-vortex condensation, is a 
local magnetic After a relaxation of $ to zero energy density the ground state 
must exhibit the associated Z 2 degeneracy. We conclude that for SU(2) the param- 
eters f = ±1 must be identified: They describe the same minimum of $'s potential. 
The parameter value f = corresponds to the other degenerate minimum. The 
center flux carried by a given flux line is associated with the differences in $'s phase 
in each minimum of <I>'s potential. 

SU(3) case: 

For SU(3) the dynamically broken symmetry is a local magnetic Z3. As a con- 
sequence, each of the three possible values f = ±1,0 describes one of the three 
possible, distinct minima of <£>'s potential. Before relaxation to zero energy density 
local jumps of $'s phase generate two distinct species of flux loops: Each associated 
with the three differences in $'s phase modulo three. (A short jump between two 
neighbouring unit roots is equivalent to a long jump into the opposite direction 
involving the third unit root as a brief stop-over.) 

4.3. The potential of the macroscopic complex scalar field 

At T Cj m dual gauge modes decouple. Moreover, a condensate of (Cooper-like) pairs 
of single center-vortex-center-antivortex loops confines fundamental electric and 
fermionic test charges. This happens because each condensed center-vortex loop 
represents an electric dipole. A condensate of such dipoles squeezes an externally 
applied electric field into a flux tube: Oppositely charged test particles are subject 
to a linear potential at large distances. Thus the center phase is confining both test 
charges and all gauge modes: There is complete confinement. (The Polyakov loop 
expectation is zero below T Cj m, the 't Hooft loop expectation <i>, which is the dual 
order parameter for confinement, jumps to a finite value.) 

The effective action for the center phase thus is only a functional of $ and $. 
Moreover, thermal equilibrium (that is, periodicity in Euclidean time) is no longer 
applicable to constrain $'s potential. According to our discussion in the last section 
3>'s potential Vc must be (i) invariant under center jumps only (invariance under a 
larger (continuous or discontinuous) symmetry is forbidden), (ii) it must allow for 
fcrmion creation by center jumps, and (iii) center-degenerate minima of zero energy 
density have to occur. Moreover, (iv) we insist on the occurrence of a single mass 
scale Ac only. From (i) we conclude that Vc can not be a function of $$ alone. 
SU(2) case: 

The only potential Vc satisfying (i),(ii), (iii), and (iv) is given by 



V c = vd~v c = -2--A c $ 




') 




') 



(173) 
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Fig. 21. The potential Vc = vq(Q)vc(&) for the center- vortex condensate <!>. Notice the regions 
of negative tangential curvature inbetween the minima. 



The zero-energy minima of Vm are at <I> = ±Ac. It is easy to show that adding or 
subtracting powers (<j> _1 ) 2 ' +1 or $ 2k+1 in vc, where k,l = 0, 1, 2, 3, • • • , generates 
additional minima and thus destroys the center degeneracy of the ground state after 
relaxation and/or violates the demand for zero energy-density at a finite value of $ 
in these minima (requirement (iii)). 

SU(3) case: 

The only potential Vc satisfying (i),(ii), (iii), and (iv) is given by 

Vc = vc-v c ^(^-^ (^-* 2 ) • (174) 

The zero-energy minima of Vc are at <I> = Acexp [±^jp] and $ = Ac- Again, 
adding or subtracting powers ($ _1 ) 3/+1 or ($) 3 ' c_1 in vc, where I — 0, 1,2,3,- •• 
and k = 1, 2, 3, 4, • • • , violates requirement (iii). 

In Fig.l2i"1 plots of the potentials in Eq. (fTT3|) and Eq. ([174]) are shown. The 
ridges of negative tangential curvature are classically forbidden: The field $ tunnels 
through these ridges, and a phase change, which is determined by an element of the 
center Z 2 (SU(2)) or Z 3 (SU(3)), occurs locally in space. This is the afore-mentioned 
generation of one unit of center flux. 



4.4. Thermodynamics close to the Hagedorn transition 

The action describing the process of relaxation of $ to one of Vm's minima is 

S = J d 4 x Q d^d^ -\v C ^j ■ (175) 

In contrast to the electric and magnetic phases the action S in Eq. (|175[) does not 
determine a classical, macroscopic ground state if $ is not in one of Vc's minima. 
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Though tunneling processes occur in real time they can be described by a Euclidean 
simulation, WKB-like approximations are thinkable. Alternatively, the computation 
of fermion creation rates can be performed on a finite-temperature lattice based on 
the theory (|175[) . An interesting object to be measured is <I>'s two-point correlator 
H(x) = ($(x)$(0)). Projecting onto a given spatial momentum p at a given temper- 
ature, the strength of intermediate tachyonic modes can be extracted by a Fourier 
analysis of the t dependence in J d 3 x exp[ip • x] n(r, x). This gives a measure for 
the density of states p n for fermions of mass ~ n Ac and spatial momentum p^l. 

Let us estimate p n for SU(2). The multiplicity of fermion states with n self- 
intersections is given by t wice the number L n of bubble diagrams with n vertices 
in a scalar A</> 4 theory. In — the minimal number of such d iagrams L n rn i n was 
estimated to be L n ^ m i n = n!3~™. Using Stirling's formula this can be approximated 
as 

L n ,min ~ (3e)~™ n n 



_ V2tt 
~ ~3~ 

3~~ 

for ii 3> 1. So the number F r 
from below roughly by 

\/2T 

f„~2x 



! -\[> n^logn — (log3 + 1) 
exp [n log n] 

of fermion states with mass 



ilogn 



exp [n log n] 



/87T 



exp [n log n] 



(176) 

nAc is bounded 
(177) 



3 r L ° ' 3 
We now estimate the density p n> Q of fermion states at rest p = (or p(E = nAc)) 
by differentiating F n with respect to n and dividing the result by the level-spacing 
8m n — Ac- We have 

/8tt 



Pnfi > 



p{E)> 



3Ac 



exp[n log n] ^log? 



1 



3Ac eXP te l0g A^ ] l l0g A^ 



1 



(178) 



Eq. (|178p tells us that the density of static fermion states is more than exponentially 
increasing with energy E. The partition function Z$ thus is estimated as 



> 



dEp(E) n F {j3E) 



/87T 



3A r 



/>oo 


' E ' 


/ dE exp 




)E* 



exp[-/3J5] 



(179) 



where E* 3> Ac is the energy where we start to trust our approximations. Thus 
Z$ diverges at some temperature Th < Ac- (Strictly speaking, Th = according 
to Eq. (|178p . This is an artefact of our assumption that all states are infinitely 
narrow. There are, however, finite widths for higher-charge states (n > 1) since 
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contact interactions exist between vortex lines and intersection points. Moreover, 
in the real world higher-charge states are even broader due to their decay and their 
mutual annihilation into charge-one and charge-zero states. This happens because 
an SU(2) theory, which is not confining at the present temperature of the Universe, 
mixes with the theory under consideration at large temperatures and thus couples 
its massless gauge mode - the photon - to the Z2 charges of the latter. The larger 
n the broader the associated state and the less reliable our assumption of infinitely 
narrow states.) This is the celebrated Hagedorn transition. (At Th the entropy 
diverges and the system condenses self-intersecting center-vortex loops into a new 
ground state: The monopole condensate of the magnetic phase. The process of 
monopole condensation from below violates the spatial homogeneity of the system: 
Z2 charges loose their mass by dense packing only. Thus the Hagedorn transition is 
genuinely nonthermal.) 

Once $ has is settled into Vc's minima $ min there are no quantum fluctuations 
<5$. Writing $ = |$| exp ij- , this is a consequence of the following fact: 



8 (su(2)) risen 

18 (SU(3)) • (180) 



Thus radial and tangential fluctuations around $ min would have a mass mj$ which 
is sizably larger than the compositeness scale |<& min | for both SU(2) and SU(3). Since 
\Ps$ + < |$min| 2 for any allowed Euclidean momentum > this means 

that the fluctuations <5<£> are absent: After relaxation the ground state of the center 
phase has a vanishing pressure and a vanishing energy density. 



5. Matching the pressure 

The mass scales Ae and Am, which determine the modulus of the adjoint Higgs 
field 4> and the moduli of the monopole condensates (p (SU(2)) and <p\, ip2 (SU(3)), 
respectively, are related. This derives from the fact that across a second-order tran- 
sition the pressure is continuous. We have 

(SU(2)), 

(SU(3)). (181) 

Across the magnetic-center transition the pressure is not continuous. We may, how- 
ever, estimate the scale Ac by assuming thermal equilibrium at the onset of this 
transition. (This assumption underlies Eq. (|170[) .) This gives 



Am = U 



A 



M 



720tt 2 
720vr 2 



1/3 



A, 



1/3 



A, 



A M ^2 1 / 3 A C , (SU(2)) and A M ~ A c , (SU(3)) . (182) 
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6. Thermodynamical quantities 
6.1. Results 

In this section we present our numerical results for one-loop temperature evolutions 
of thermodynamical quantities throughout the electric and magnetic phase. 

SU(2) case: 

In the electric phase the ratio of pressure P and T 4 is given as 



P 



A| 



(2tt)< 



(2P(0) + 6P(2a)) +2AS 



(183) 



where the function P(a) and the dimensionless mass parameter a are defined in 
Eq. (|100p and Eq. (j9"9")l , respectively. In the magnetic phase we have 



P 



(27r)6 naj+ 2 . 



(184) 



where a is defined in Eq. (I149[) . 

In the electric phase the ratio of energy density p and T 4 is given as 



P_ 



2A 4 

--^(2p(0) + 6p(2a)) + 2A i 



where the function p(a) is defined as 



dx x' 



\/x 2 + , 



exp(v'x 2 + a 2 ) — 1 



In the magnetic phase we have 

(2^) 4 



_P_ 

rp4 



\ 4 



6A M -/ s , Am 



The ratio of entropy density and T 3 is given as 



s \ . 
j^ = jqip + P) 



(185) 



(186) 



(187) 



(188) 



Because the ground-state contributions cancel in 7^3 this quantity is not as infrared 
sensitive as, e. g., ^or^: Lattice simulations are in a position to correctly measure 
the entropy density at low temperatures. 

SU(3) case: 

In the electric phase we have 



P 



and 



P_ 



i^r 



A| 



2A| 

(2tt)« 



2A| 
(2tt)« 



(4P(0) + 3(4 P(o) + 2 P(2o))) + 2Ai 



(4p(0) + 3(4 p(o) + 2p(2a))) + 2A; 



(189) 



(190) 
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Fig. 23. y% as a function of temperature for SU(3) as obtained on a (16 3 X 4)-lattice using the 

differential method wit h a universal two-loop pe rtur bative f3 function EBE51 anc j using the integral 
method (solid line) ESI. The fi gure is taken from 12^1 



In the magnetic phase we have 

P _ (2tt) 4 



and 



ji4 



P_ 

ji4 



\ 4 



(191) 
(192) 

jggHD where now the SU(3)- 



The ratio of entropy density and T 3 is given in Eq. 
expressions for P and p have to be used. 

The result for ^ is plotted in Fig. [22] as a function of temperature throughout 
the electric and magnetic phase, Fig. |2"31 depicts SU(3)-lattice results. Notice that 
the pressure is negative in the electric phase close to Xe.c and even more so in the 
magnetic phase where the ground state strongly dominates the thermodynamics of 
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infrared sensitive quantities. Notice also the negative pressure in Fig. [23] obtained 
close to the phase transition when the differential method is used in the lattice 
simulation. (For a discussion of differential versus integral method see Sec. 16. 2. 2] ) 
We conclude that the finite-size constraints of realistic lattices have a severe effect 
on the obtained values of the pressure shortly above the first confining transition 
and even more so below this transition. 

Let us now discuss the behavior close to T Ct M where the thermodynamical rela- 
tion 

dP = SdT (193) 

starts to be violated. Eq. (|193|) implies that in a homogeneous, thermalized system 
the pressure needs to be a strictly monotonic increasing function of temperature 
since S is never negative. Crossing the point T c m from above, the pressure jumps 
from a negative to a positive value. (On the magnetic side of the phase boundary 
the ground state strongly dominates the excitations, on the center side the vortex- 
condensate has zero pressure while fermionic excitations give a positive contribu- 
tion.) There are two ways of seeing that thermal equilibrium must break down close 
to T Ci jvf. First, spatial homogeneity starts to be badly violated by discontinuous and 
local phase changes of the field $ as soon as the system starts to condense center- 
vortex loops. The derivation of Eq. (|193[) . however, relies on thermal equilibrium 
and thus on spatial homogeneity. Second, one may assume thermal equilibrium and 
then lead this assumption to a contradiction. In thermal equilibrium the spectral 
density p(E) in the center phase is more than exponentially increasing with energy 
E, see Eq. (|178[) . Thus the partition function diverges at T = Th < Ap: A homo- 
geneous system would need an infinite amount of energy per volume to increase its 
temperature beyond Th- But this is a contradiction to the fact that a magnetic 
phase exists for T > Th- (In an extended thermalized system the transition from 
the center phase to the magnetic phase is accomplished by an increase of the over- 
all energy density: The excitation of very massive dual gauge modes, though very 
unlikely, is furnished energetically by the large energy residing in the system. If 
the considered system is of a small spatial extent, such as the interaction vertex in 
a scattering process, then the total energy of the system, e.g., the center-of-mass 
energy being deposited into a vertex, needs to be larger than the very large mass of 
the dual gauge modes on the magnetic side of the phase boundary.) Thus thermal 
equilibrium breaks down for T ~ Th- 

The result for ^ as a function of temperature throughout the electric and mag- 
netic phase is shown in Fig.[M] Fig. [25] depicts an SU(2)-lattice result Notice 
the (small) discontinuities at \ c ,e- Their occurrence is explained by the fact that 
by crossing the electric-magnetic phase boundary the system needs to generate an 
extra polarization for each dual gauge mode compared to the two polarizations of 
a TLM mode. Extra polarizations are extra fluctuating degrees of freedom which 
increase the energy density on the magnetic side of the phase boundary. The situa- 
tion is somewhat peculiar: On the one hand, the order parameter for the dynamical 
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Fig. 24. ^ as a function of temperature for SU(2) (left panel) and SU(3) (right panel). The 
horizontal lines indicate the respective asymptotic values, the dashed vertical lines are the phase 
boundaries. 

breaking of the dual gauge groups U(1)d (SU(2)) and U(l)|, (SU(3)) is contin- 
uous. On the other hand, there is a small amount of latent heat being released 
across the magnetic-electric transition. (That is, by heating up the system starting 
in the magnetic phase.) As we shall see in Sec. [3 this is the reason for a dynam- 
ical stabilization of the temperature of the cosmic microwave background against 
gravitational expansion. (Thus we may look forward to enjoy the privilege of the 
photon's masslessness for another sizable fraction of today's age of the Universe 
^U) Again, the energy density is dominated by the ground-state contribution in 
the electric phase close to the electric-magnetic transition and even more so in the 
magnetic phase. Notice also that p = —P at the point, where the system starts to 
condense center-vortex loops, that the magnetic phase is narrower for SU(3) than 
it is for SU(2), and that ^ dips in a much steeper way at the electric-magnetic 
transition for SU(3) than for SU(2). 

The result for the interaction measure ^ = is shown in Fig. [521 Figs.l2"Tl 

and [28] are lattice results. Notice the rapid approach to the free-gas limit in Fig. [26] 
and the large values of ^ in the magnetic phase. Interestingly, there is a small 
bump to the left of the phase boundary in Fig. [27] 

The result for the ratio of the specific heat per unit volume cy = 4& and T 3 is 
shown in Fig.l2l?l Fig.lgQlis an SU(2)-lattice result^. The quantity peaks both 
at the electric-magnetic and the magnetic-center transition. The finite peak at the 
former phase boundary is in agreement with the electric-magnetic transition being 
second-order. Moreover, we have SSL ~ 3 §£■ L, . This explains why 

i IT C-E ;su(3) I J IT CiE ;su(2) 1 ■' 

lattice simula tions p refer to identify the confining transition as weakly first order 
for SU(3), see anc i references therein. Now, 3 7^ 00 but in the vicinity of T c e 
lattice results for infrared sensitive quantities, such not reliable anyway. 
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Fig. II, The ratio e/e SB for a 10 3 X 3 lattice versus temperature. The full points are icosahedral, the 
crosses full SU(2) group Monte Carlo results; the line is a spline fit to the data. 

Fig. 25. as obtained from the SU(2)-lattice simulation in LZSJ. 
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Fig. 26. The interaction measure ^ as a function of temperature for SU(2) (left panel) and 
SU(3) (right panel). The asymptotic value in both cases is = 0, the dashed vertical lines are 
the phase boundaries. 



The result for ^ is shown in Fig. [22 Fig. [321 depicts a lattice result for SU(3) 
obtained with the differential method^!. The entropy density S is a measure for the 
'mobility' of gauge modes. Notice the jump of S/T 3 which, again, is explained by the 
additional polarization of the dual gauge mode in the magnetic phase. Notice also 
that tpr vanishes at the point T c m where the system condenses center vortices. At 
this point dual gauge modes are infinitely heavy: The thermodynamics is completely 
determined by the ground state. The numerical agreement between the lattice result 
(d) (largest lattice) in Fig. [32] and the SU(3)-result in Fig.[3T] is striking. The two 
data points to the left of the jump in (d) indicate that an ambiguity exists for the 
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Fig. 15. The difference (e — "V)/7"' versus temperature for a 10 3 X 3 lattice. 
Fig. 27. as obtained from the SU(2)-lattice simulation in ESI 

value of ;Jj very close to the transition. The jump itself corresponds to the large 
slope of on the electric side of the phase boundary in Fig.[3TJ The observed 
agreement is explained by the small sensitivity of the quantity ^ on the ground- 
state physics making the finite- volume lattice simulation reliable close to the electric 
magnetic transition. 

The data files needed to generate the plots in Figs. OH EMI EH US and [SH are 
provided by the author upon request. 



6.2. Comparison with the lattice 
6.2.1. Specific observations 
SU(2) case: 

The results of an early lattice measurements of the energy density p and the interac- 
tion measure A = p — 'SP in a pure SU(2) gauge theory were reported in . In that 
work the critical temperature T c for the deconfinement transition was determined 
from the critical behavior of the Polyakov-loop expectation and the peak position of 
the specific heat using a Wilson action. The function A(T) was extracted by multi- 
plying the lattice (3 function with the difference of plaquette expectations at finite 
and zero temperature ( sym metric lattice). This assures that A vanishes for T — > 0. 
What is subtracted in at finite T is, however, not the value A(T = 0) since 
the plaquette expectation on the symmetric lattice is multiplied with the Gnite-T 
value of the (3 function. Apart from this approximation, the use of a perturbative 
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FIG. 2. (a) The pressure versus T /T c for /V,. = 4, 6, and 8, 
integrating the interpolations for the action density. Solid 
curves show our parametrization. The dashed curve for N ^ = 4 
is the result of using the temperature scale from the effective 
coupling scheme. The horizontal dash-dotted lines show the 
ideal gas values for JV^ — 4, 6, and 8; the horizontal dashed 
line is the continuum value. In (o) we show the difference 

Fig. 28. ^3- and t^j- as obtained from the SU(3)-lattice simulation in 



(3 function was assumed for all temperatures. The simulation was carried out on a 
(rather small) (10 3 x 3)-lattice. 

Let us compare our results with those of ^1 The lattice results for p in ^3 
differ drastically from our results for temperatures close the first confinement, that 
is, the electric-magnetic transition. (The lattice is doomed to produce incorrect 
results for infrared sensitive quantities close to the electric-magnetic transition and 
in the magnetic phase: A finite spatial lattice size L cuts off physical correlations on 
length scales > L since the correlation length Im — Mm 1 > L close to the electric- 
magnetic transition and Im — oo in the magnetic phase. Here M m denotes the mass 
of a magnetic monopole.) 
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Fig. 29. ^ as a function of temperature for SU(2) (left panel) and SU(3) (right panel). The 
horizontal lines signal the respective asymptotic values, the dashed vertical lines are the phase 
boundaries. 
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Fig. 12. The specific heat per unit volume divided by T 3 versus temperature, obtained from the fit in fig. 

1 1 to £/e SB . 



Fig. 30. Tpj as obtained from the SU(2)-lattice simulation in 



75 



We obtain 



P 

PSB 



1.27, 



(194) 



T~1.5T C 



2 . 

where p SB = ^-T denotes the Stefan-Boltzmann limit (ideal gas of three species of 
massless gluons with two polarizations each). On the lattice this ratio is measured 



to be smaller than unity: 

■' PSB 



T~1.5T C , E 
P_ 

T~5T C , E 



PSB 



= 0.84. At T ~ 5T C ^ E we obtain 
~ 1.33 (195) 
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(a) (b) (c) (d) 

Fig. 32. 7^3 as a function of /3 obtained in SU(3) lattice gauge theory using the differential method 

and a perturbative beta functionl^§l The simulations were performed on (a) 16 3 X 4, (b) (24 3 X 4)-, 
(c) (16 3 X 6)- (open circles) and (20 3 X 6)- (closed circles), and (d) (24 3 X 6)-lattices. Using the 
(24 3 X 6)-lattice, the critical value of is between 5.8875 and 5.90. 



while the lattice measures — = 0.85. 

PSB T~5T C , E 

Our results for the pressure P are negative for T close to T CiE (see Fig 
much in contrast to the positive values obtained in At T ~ 5T c e we obtain 



P 

PsB 



1.31 



(196) 



T~5T CjE 



while the lattice measures 



~ 0.88. (On the lattice P is extracted from 

T~5T CjE 

the measured values of A and p, and P SB = j^T 4 denotes the Stefan-Boltzmann 
limit.) 

Notice that the results in Eq. (|195p and Eq. (|196|) are very close to the ratio 
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2X g +2 1 of the number of degrees of freedom in a gas of two species of (nearly) 
massless gluons (three polarizations per species) and one massless species and a gas 
of massless gluons. (At As ~ 75 the value of the mass parameter is a ~ 27r|g| ~ 
0.086. Thus the Boltzmann suppression is small for TLH modes, compare also with 
Fig.1261) At extremely high temperatures a TLH mode 'remembers' its massiveness 
at low temperatures in terms of an extra polarization. The latter originates from a 
tin y m ass which solves the infrared problem of loop expansions, for formal arguments 

The peak- value of the specific heat is about §&| T E ~ 20 while it is measured 
to be ~ 8 on the lattice. Moreover, we have y?| T ~ 4.8 while the lattice obtains 
a value ~ 0.85. The much lower values obtained on the lattice are not surprising: 
Finite lattice sizes cut off existing long-range correlations at TV 

No result for the entropy density was directly reported in ^ 

SU(3) case: 

Here we discuss the results obtained in with a Wilson acti on on the lattice of 
the largest time extension, Np — 8, and the results obtained in l ^9 | 78 | 

In the vicinity of the transition point T Ci e the situation for both p and P is 
similar to the SU(2) case: Drastic differences between the lattice measurements and 
our results occur. Again, P is negative close to T Cj £ contradicting the positive val- 
ues obtained with the integral method in A lattice simulation ^1 of P, which 
used the differential method, has reported negative pressure for T shortly above the 
transition already in 1988. The most negative value of P/T 4 ~ —0.5 obtained in^l 
very close to the phase transition is down by a factor of about 0.19 as compared 
to our result at the electric-magnetic transition, see Figs.l22l andl23l Again, this is 
explained by the finite lattice-size cutoff on ph ysic al long-range correlations. The 
lattice- result obtained with the integral method LLLil is by construction positive defi- 
nite, see Sec. l6.2.2| and thus it is unphysical. For that reason we renounce a (useless) 
comparison of our res ult s for pressure, energy density, and interaction measure with 
those obtained in ^1. In ^1 nly the dependence of ^ on the lattice coupling was 
presented. One can use Fig. [23] and the temperature dependence of in Fig. |2"51 to 
auge part icular values of this quantity against temperature. (In both simulations 
and ^1 t ne universal part of the two-loop perturbative (3 function was used to 
relate lattice coupling to lattice spacing.) For example, a value of £^ ~ 1.6 in^l 
corresponds to a value A ~ 1.5 inlZSl The latter is associated with a temperature 
T = 3.2T c ,e by virtue of Fig.El We have 

P 



-2.2. (197) 

3.2T C , E 



This is larger than the result obtained in and explained by the insufficient account 
of infrared correlations in the lattice simulation. These correlations generate masses 
for six out of eight gluon species, thus extra polarizations, and therefore a larger 
value for 
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Our asymptotic (Xe = 35) values for P and p are 



P 



where psB = 



15 



P S B 
7T 2 T 4 = 



1.30. 



P 

PSB 



1.37. 



(198) 



3Psb- Both values in Eq. (I198|) are close to the ratio 



R = ^ = 1.375 of the numbers of polarization in a free gluon gas, where six gluon 
species have a tiny mass, and in a free gluon gas where all gluon species are massless. 

The entropy density approaches zero for T \ T Cj m, see Fig. ([31]) . This expresses 
the fact that dual gauge modes are decoupled (infinite masses): The ground state 
strongly dominates the thermodynamics. 



6.2.2. Differential versus integral method 

What are t he rea sons for the qualitative difference betwee n the pressure-results 
obtained in l ?6 | 77 | using the integral method and in using the differential 

method? While the differential method is based on the definition 

which is proper for a lattice of finite volume V, the integral method assumes the 
thcrmodynamical limit V — > oo from the start. In this limit one has 

P = , (200) 

and thus the pressure equals minus the free energy density. In Eqs. (|199p and (|200[) 
Z denotes the partition function. 

The official reason for the introduction of the integral method, see for example 
EH! W as that one wanted to avoid the use of the imprecisely known (3 function 
in the strong-coupling regime of the theory. (Based on the definition in Eq. (|199[) , 
the (3 function multiplies the sum of spatial and time plaquette averages in the 
expression for the pressure.) When using the definition in Eq. (|200[) . the derivative of 
the pressure with respect to the bare coupling (3 {(3 = % for SU(3)) can be expressed 
as an expectation over minus the sum of spatial and time-like plaquettes without 
the beta-function prefactor. Thus the pressure is, up to an unknown integration 
constant, obtained in terms of an integral of a sum of plaquette averages over (3. 
The integration constant is chosen in such a way that the pressure vanishes at a 
temperature well below T c . Instead of only integrating over minus the sum of spatial 
and t ime-like plaquette expectations an extra term was added to the integrand 
| 76 | 77 | to assure that the pressure vanishes at T — 0. The added term equals twice 
the plaquette expectation taken on a symmetric lattice (the expectation at T = 
). We stress that this prescription does not follow from the definition in Eq. (|200|) . 
Moreover, the assumption that P = for T ~ 0.8 T c or so is a strong bias. (There 
are massless fermionic particles in the center phase which keep the total pressure 
positive for temperatures comparable to this value.) 
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The results for P(T) obtained when using the integral method show a rather 

large dependence on the spatial size and the time extent N T of the lattice ^51. We 

believe that this reflects the considerable deviation from the assumed thermody- 

[771 

namical limit for realistic lattice sizes. The problem was addressed in L^-*- 1 where a 
correction factor r was introduced to relate P, obtained with the integral method, 
to P, obtained with the differential method. For a given value of N T the factor r was 
determined from the pressure-ratio at g — 0. Subsequently, this value of r was used 
at finite coupling g to extract the spatial anisotropy coefficient c CT (essentially the j3 
function) by demanding the equality of the pressure obtained with the integral and 
the differential method. In doing so, twice the plaquette expectation at T — was, 
again, added to minus the sum of spatial and time-like plaquette expectations in 
the differential-method expression for the pressure. It may be questioned whether 
a simple correction factor r does correctly account for finite-size effects and, if yes, 
whether it is justified to determine r in the limit of noninteracting gluons. (The ca- 
valries obtained in this way do not coincide with those obtained in^H) In addition, 
it seems that the imprecise knowledge of the (3 function, which contains informa- 
tion about fluctuations in the ultraviolet, is much less of a problem for a lattice 
simulation of the pressure than the missing infrared physics is (finite lattice size) 

M 

Using the universal part of the two-loop perturbative beta function in the dif- 
ferential method, negative values for the pressure were obtained for T close to T c 
j n [79| Moreover, a rapid approach of p and P to their respective free-gas limits was 
observed. This is in qualitative (but not quantitative) agreement with our results, 
see Figs.Ha [H and [22 

7. Implications for particle physics and cosmology 

In this section we provide outlooks on the implications of the nonperturbative ap- 
proach to SU(2) and SU(3) Yang-Mills thermodynamics in view of so-far unex- 
plained phenomena in particle physics and cosmology. The way of how selected 
problems are addressed in this section is preliminary, mostly qualitative and thus 
should not be understood as the final say on the matter. Rather, we try to provide 
a certain amount of stimulus for future developments. 

7.1. A Planck-scale axion: Cosmic coincidence today and CP 
violation 

Among the gauge groups SU(iV) (N finite) we regard SU(2) and SU(3) as particular 
due to their unique phase diagrams. We have come to appreciate that nature seems 
to prefer situations with a unique outcome. Thus we tend to believe that dynamics 
subject to a finite gauge symmetry, that is, dynamics below the Planck scale Mp ~ 
1.2 x 10 19 GeV, obeys the SU(2) or SU(3) gauge principle. A possible scenario would 
be that at Mp an SU(iV) gauge symmetry (N = oo) is dynamically broken into a 
four-dimensional low-energy manifestation involving several SU(2) and SU(3) factor 
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groups, which can behave in an electric-magnetically dual way to one another, and 
into nonfluctuating gravity. 

This set-up may, in fact, be described by the low-energy sector of a bosonic string 
theory whose vacuum inst ability is resolved in terms o f tachy on condensation in the 
presence of a D-brane^l A low-energy Kaluza-Klein IMI65| compactification of the 
Weyl-invariant and thus (i-dimensional bosonic string theory (in a flat background 
d = 26) to four dimensions yields gauge symmetries which are associated with the 
isometries of the d— 4 dimensional compactification manifold: We assume that these 
isomctries are products of SU(2) and SU(3) c orresp onding to a compactification 
manifold which, locally, is S3 x (53 x S5) x • • • I 60 | 59 [ j n addition, there are a low- 
mass scalar dilaton field ip and a massless antisymmetric tensor B^. While the 
former may drive the monopole condensation process within, say, an SU(2) theory 
of Yang-Mills scale Ai ~ Mp and thus may trigger the inflation of the Universe 
(making it spatially fl at) the latter may be responsible for adiabatically generated 
density perturbations i2 - Li . 

The scales of the SU(2) and SU(3) factors would dynamically bet set into a 
certain hierarchy: Ai ~ Mp > A2 > ■•• > A CMB . The scale A CMB is associated 
with an SU(2) theory that is not confining at the present temperature of the Uni- 
verse: Being at the electric-magnetic transition this theory generates the photon as 
its only massless excitation, see Sec. 17.21 Being in its center phase at a temperature 
T ~ Ai ~ Mp, SU(2)i generates fermions by re-heating (single and self-intersecting 
center-vortex loops). Because SU{2)\ used to be part of SU(A = 00) and thus was 
mixing with the other factors these fermions couple to the gauge fields of these fac- 
tors. Since higher-charge states (self-intersecting center- vortex loops) with opposite 
charges are generated in equal amounts we would expect that they quickly anni- 
hilate into charge-zero states (no self- intersections). (Re-creation of higher-charge 
states after annihilation is less likely due to the redshift of the spectrum by the rapid 
power-law expansion of the Universe furnished by a free-gas equation of state. The 
latter originates from the gauge- mode excitations of the other factors). Although 
the massless fermions (single center-vortex loops) do not couple to the propagating 
gauge fields of the other factors by naked gauge charges they do so by their dipolc 
moments. Considering these massless fermions to be fundamental, there is a global, 
axial U(l) symmetry which, howev er, is d ynamically broken I^Hl and anomalous due 
to the calorons of the other factors I 8Q | 90 |_ 

Integrating out SU(2)i's massless fermions, the relevant composite field is a 
(canonically normalized) axion field a whose coupling to the gauge fields of the 
other factors is 

CMB 

£*u.a. = ^ ^ tr F^t/. jF^u.i (201) 

where F ~ Alp denotes the Peccei-Quinn scale The scale at which SU(2)x'b 
fermions come into existence and which measure the magnitude of the Cooper-pair 
condensate involving the massless species. The sum is over the other factors, and 
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Fixu,i — 5 F a fj :i is the dual field strength. Eq. (|201[) represents a term in the 

action for the other factors (deconfining at T ~ Mp) which violates parity (P) and 
charge conjugation (C) symmetries. 

While a would be a massless phase if the axial anomaly was absent and the axial 
U(l) only was broken dynamically the anomaly-induced coupling in Eq. (|201[) gives 
rise to an axion potential V a = Y^SS Va,i- The operator trF M „ ^F^i measures 
the average topological charge density carried by propagating gauge fields. This 
is a conserved quantity for T > A„ and thus it is independent of temperature. 
Integrating over topologies and accounting for dimensional transmutation we have 

f (l — cos j?) Af , (theory i in electric or magnetic phase) , (202) 
a ' 1 \ (theory i in center phase) . 

For a smaller (but not much smaller) than F the cosine in Eq. (|202| can be expanded, 
and the axion mass-squared at T ~ Mp reads 

CMB CMB i 4 

i=2 i=2 p 

As the temperature of the Universe falls below A2 the associated theory fails to 
contribute to the axion mass-squared and so forth. 

Now in an expanding Friedmann-Robertson- Walker Universe the (spatially ho- 
mogeneous) axion field a satisfies the following equation of motion 

H + 3Ha + mla = 0. (204) 



The Hubble parameter H = ^ = y 8 ^[P is determined by the energy density 

P(F) = ^^=2 Pi(F) of the Universe. At T ~ Mp this energy density is given by the 
Stefan-Boltzmann limit of the theories i = 2, • • • , CMB if the hierarchy between Mp 
and A 2 is sufficiently large. (The SU(2) theory with i = 1 went center, subsequently 
experienced a strong dilution of its massive excitations, and does not produce a 
ground-state contribution to p at T ~ M P .) We have p l (T) = ^-T A (SU(2)) 
and Pi (T) = T 4 (SU(3)). Thus H is radiation-dominated at T - Mp and 
H 2 ^> m 2 . But this means that a is frozen to the slope of its potential with an 

amplitude a ~ Mp. For T ~ A 2 the Universe's energy density remains radiation- 

A 4 

dominated. The axion mass-squared, however, is reduced by the value j^- since the 

theory with i = 2 went center. If A 2 <C Mp then the mass-squared ~ Y>i=3 °^ 
the axion generated by the theory with i = 2 is much larger than H 2 . By Eq. (|204p 
this means that this axion rapidly relaxes to the minimum of its potential and thus is 
irrelevant for subsequent cosmology. Moreover, the fermions generated by the center 
transition of the theory with i = 2 exhibit a large asymmetry in fermion number. 
This is qualitatively true since all three Sakharov conditions are satisfied: (i) the 
center transition is nonthermal (Hagedorn), (ii) there is a local violation of fermion 
number since fermions are nonlocal objects, and (iii) the generation of fermions 
takes place in the presence of CP violation (the frozen- in Planck-scale axion a). 
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Fig. 33. The fate of a Planck-scale axion along the Universe's evolution. At temperatures sizably 
larger than Tcmb the axion is frozen to the slope of its potential by cosmological friction (H m a ), 
for T ~ Tcmb axion mass and Hubble parameter become comparable: The axion starts to roll down 
its potential. 

This goes on until T — Tcmb where the last theory in the chain, SU(2)cmb? 
is close to its center transition (more specifically, at the electric-magnetic phase 
boundary today), see Fig. [33] Here H 2 is dominated by the ground-state energy. 
We have 



In Eq. (|205p pv denotes the energy density associated with the value of the axion 
potential at T CMB , and px is an energy density due to axion rolling. Both contribu- 
tions are comparable since H is not much larger than m a , compare with Eq. (|204p . 
The formerly frozen-in axion field a slowly starts to roll down its potential. While 
pv has an equation of state pv — —Pv the kinetic contribution px is associated 
with an equation of state Pk = which is the same as that of nonrelativistic matter. 

Cosmic coincidence may have an explanation in terms of a Planck-scale axion 
and an SU(2) Yang-Mills theory of scale A nMR comparable to the present temper- 
ature of the CMB. For related ideas see ^21 The alert reader may object that the 
9 angle in Quantum Chromodynamics (QCD) is constrained to be an extremely 
small number by a measurement of the neutron's electric dipole moment. On the 
other hand, the mass of the r{ is much larger than the pion mass. Thus one would 
conclude that the Planck-scale axion is irrelevant in the former while it is relevant 
in the latter case. What is the resolution of this puzzle? The electric dipole moment 
of the neutron is measured with a photon of momentum much smaller than the 
QCD confinement scale. Thus this photon does not probe a phase of QCD with 
propagating gauge bosons: The operator F^F^ has a vanishing expectation. The 
rf , on the other hand, is generated in a scattering process involving propagating 
gluons: Inside the vertex the operator F^F^y has a finite expectation. 

We will see in Sec. 17.21 that the ground-state contribution of SU(2) CMB in the 
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"CMB 
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absence of the Planck-scale axion is small in comparison with the measured value of 
today's cosmological constant. The scenario outlined above does not yet explain the 
origin of clustering dark matter as it is observed in the anomalous rotation curves 
of galaxies but we will see below that the decoupled TLH modes of SU(2) CMB are 
candidates for this form of matter. 

7.2. Electroweak sector of the Standard Model: Nature of leptons, 
electroweak symmetry breaking, masses of intermediary vector 
bosons, intergalactic magnetic fields, and solar wind 

Here we would like to propose a formulation of the electroweak sector of the Stan- 
dard Model in terms of pure SU(2) Yang-Mills theories. 

SU(2) C MB • 

Let us first discuss the U(l)y factor of the electroweak gauge group SU(2) w xU(l)y. 
We claim that this factor is the unbroken subgroup of an SU(2) Yang-Mills theory 
of scale comparable to that of the CMB temperature: T CMB = 2.728 A" = 2.351 x 
10 -4 eV. Only one point in the phase diagram of this theory exists, the boundary 
between the electric and magnetic phases, where this claim is in accord with obser- 
vations: The photon is unscreened and practically massless (m 7 < 10~ 14 eV from a 
precision measurement of the Coulomb potential see for a discussion on why 
stronger bounds are unreliable), see Fig. [HI Thus we identify T CMB = T c _e- Notice 
that isolated charges in the electric phase of SU(2) CMB have a dual interpretation: 
What is a magnetic monopole in SU(2) is an electrically charged particle w.r.t. 

U(l)y. 

The energy density p gs of the ground-state at T CMB is p gs — 4:TtT cmb A 3 e = 
4tt x (2.351 x 10~ 4 eV) A| = 2A c ,_eA| = 27.7 A|. Moreover, we have A E = 1.065 x 
10~ 4 eV. Substituting this into p gs , we have p gs = (2.444 x 10~ 4 eV) 4 . This is about 
0.36% of the commonly accepted value of today's dark energy density (10 _3 eV) 4 . 
The dominating, missing part would be generated by a slowly rolling Planck-scale 
axion, see Sec. 17. II 

An immediate question to answer is why the masslessness and the unscreened 
propagation of the photon is a singled-out situation. The answer to this question is 
encoded in Fig.1241 The energy density of an SU(2) Yang-Mills theory dips at the 
electric-magnetic phase boundary. On the electric side this is explained by the ther- 
modynamical decoupling of TLH modes, on the magnetic side an extra polarization, 
which costs energy, needs to be generated for the photon. To facilitate the jump 
in energy density for 8U(2)cmb to reach the magnetic phase thermal equilibrium 
needs to be violated by the eventually fast rolling Planck-scale axion. The photon 
acquires mass and the ground state becomes superconducting (electrically charged 
monopoles condense) . It is suggestive that the occurrence of intergalactic magnetic 
fields is related to the Universe being slightly out of thermal equilibrium due to 
(slow) axion rolling today. 

At Tqmb TLH modes decouple thermodynamically. Recall that their mass is 
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given by m TLH = 2e|0| and e^ec = oo at T c> e = T CMB . In the real world e^ec is large 
but not infinite because SU(2) CMB is not the only Yang-Mills theory in the Universe. 
Taking e^ec ~ 10 6 , say, the mass of a decoupled TLH mode is m TLH ~ 57 eV. (The 
reason why we chose this value for edec is mildly justified by our discussion of the 
gauge group SU(2) e below.) Since the two TLH modes do not interact and thus are 
stable (no decay into light fermions is possible because SU(2) CMB is not in its center 
phase yet) they yield a tiny contribution to clustering dark matter. 

Finally, we would like to make a remark concerning the observed large-angle 
anomaly in the temperature-(electric)polarization cross correlation seen by WMAP 
ESI. The standard explanation is that this effect is generated by CMB photons 
scattering off electric charges which are released by an early re-ionization of the in- 
terstellar medium at redshift z ~ 10 — 20. We would like to propose that CMB pho- 
tons scatter off electrically charged and dilute monopoles at temperatures > T CMB 
which are condensed at T CMB . Since static magnetic (electric with respect to SU(2)) 
fields are completely screened in the photon propagator deep inside the electric 
phase of SU(2) cmb, see Eq. (|112j) . the effect should be weaker in the temperature- 
(magnetic)polarization cross correlation. 

SU(2) e x SU(2) M x SU(2) T : 

To relate the existence and the interactions of leptons, as they are described by the 
electroweak sector of the Standard Model (SM), of which Quantum Electrodynam- 
ics (QED) is an integral part, with pure SU(2) Yang-Mills dynamics is motivated 
by the following observations, (i) The masses of charged leptons are unexplained 
parameters in the SM. In particular, their small values on the scale of their ap- 
parent pointlikcncss is unexplained, (ii) No deeper explanation for the value of the 
magnetic dipole moment of a charged lepton other than that following from the 
Dirac equation and small radiative corrections is given, (iii) There are experimen- 
tal indications that scattering processes involving the electron or the positron do 
not obey the QED predictions if the momentum transfer is close to the mass of 
a charged lepton. (iv) No Higgs particle has been observed experimentally up to a 
hypothetical Higgs mass of ~ 115 GeV suggesting that electroweak symmetry break- 
ing takes place by a different mechanism than assumed in the SM. (v) The naive 
ground state of the SM generates a cosmological constant which is many orders of 
magnitude larger than the observed value. 

Points (i), (ii), (iv), and (v) are undisputed facts. To see that there is some 
truth to point (iii) we present experimental results. In Fig. [34] a plot of the ratio of 
experiment to theory of the wide-angle e + e~ pair-creation cross section through 7 
scattering off of the field of a carbon nucleus is shown as a function of the invariant 
mass M of the created lepton pair. Notice the substantial deviation from unity for 
M ~ 2 m M • • • 4 m M ~ (210 • • • 420) MeV. Notice also that for M > 500 MeV theory 
and experiment do agree. A much more drastic deviation within the same kinematic 
regime was seen earlier in^! This, however, was not confirmed in^l 

In Fig. [35] a plot of the differential cross section for e~e~ scattering (electron in- 
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Fig. 34. Ratio of experiment to theory of the wide-angle e + e — pair-creation cross section by 7 
scattering off the field of a carbon nucleus as a function of the invariant mass M of the lepton 
pair. Plot taken from 12ZI 




FIG. 5. Comparison with theory of the measured absolute Flo. 6. Comparison with theory of the measured absolute 

differential electron-electron scattering cross section at n=0.5 as differential electron-electron scattering cross section at w=O.S as a 

a function of the kinetic energy of the incident electron taken with function of tile kinetic energy of the incident electron taken with 

the 270° apparatus. the ISO" apparatus. 

Fig. 35. M0ller scattering of electrons, for an explanation see text. Taken from 1221 

cident on an atomic electron of a target) at a fixed fraction v — 0.5 of the incident ki- 
netic energy E kin transferred in the collision is shown for 0.6 MeV < E kin < 1.7 McV 
taken with a 270° apparatus (left panel) and for 0.6 MeV < E kin < 1.2 MeV taken 
with a 180° apparatus (right panel) The solid line indicates the theoretical re- 
sult obtained by using the M0ller formula. Notice the agre emen t with QED for large 
values of i? kin (in particular in the left panel of Fig.l3"5j). In 020] the differential cross 
section for M0ller scattering was measured for _E kin = 15.7 MeV as a function of the 
scattering angle and found to be in agreement with QED on the 0.4% error level. 
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The disagreement at the lowest value E kin — 0.6 MeV in Fig. [35] which corresponds 
to a center-of-mass energy of about 2.5 m e , is conspicuous. 

If the Yang-Mills scales of the factors SU(2) e , SU(2) M , and SU(2) T are about 
m e , m M , and m T , respectively, then the masses of the charge-one states in the 
center phase of each theory, see Fig. [20] are determined to be these values. Since 
T CMB <C m e , m^m-r these theories are in their center phases. By looking at Fig. [20] 
a g-factor of two is imperative by the asignment of angular momentum one half 
in the presence of one unit of (electric) center flux in the vortex loop. The latter 
generates the lowest nonvanishing quantum of magnetic moment (Bohr magneton) . 
Notice that Fig. [20] provides for an intuitive manifestation of the concept of spin- 
1/2: Inside the intersection core the center flux is diverted to the right above and to 
the left below such that an eddy is generated. The latter carries the electric charge 
of the soliton. 

Let us now give some qualitative arguments why the gauge group 
SU(2) CMB xSU(2) e xSU(2) M xSU(2) T together with a Planck-scale axion may be a 
viable candidate to describe the phenomenology of electroweak interactions. Why 
do the photon (SU(2) CMB ) and the massive intermediate vector bosons (magnetic 
and electric phase of SU(2) e ) couple to the charge and/or the magnetic moment 
of charge-one and charge-zero states? The answer to this question is rooted in the 
symmetry breakdown at T ~ Mp where the factors SU(2) CMB , SU(2) e , SU(2) M , and 
SU(2) T were generated out of one large gauge group whose gauge bosons where 
mixing. (In the SM the mixing of the 'photon' of U(l)y with that of SU(2)w is 
parametrized by the Weinberg angle 9 with sin 2 9 — 0.23.) The interaction of the 
photon of SU(2) cmb with the electrically charged soliton of, say, SU(2) e thus is 
furnished by an adiabatic rotation into the (massive) photon of SU(2) e when ap- 
proaching the charge of the latter and an adiabatic back-rotation into the photon of 
SU(2) CMB after the interaction has taken place. Where are the higher charge states? 
These states are instable by repulsion mediated by the photon of SU(2) CMB , and 
thus they are very broad. The density of these states, however, is over-exponcntially 
rising. Why do we only see the structure of a charged lepton in scattering experi- 
ments with a center-of-mass energy y/s comparable to the mass of the lepton, see 
Figs. 1341351 ? Radial excitations of a 't Hooft monopole have been investigated in 
The first excited level is comparable to twice the mass of the monopole ground state. 
This must semi-quantitatively also hold for a Z^ monopole (self-intersection of a 
center- vortex loop). For y/s -C rn e , m M , m T the Z2 monopole is not excitable, a QED 
point-particle description holds. For mz S> y/s S> m e , m M , m T , where mz ~ 90 GeV 
is the mass of the Z boson, the energy deposited into the vertex is converted into a 
large entropy carried by the Hagedorn spectrum of instable states. The latter protect 
the Z2 monopole against radial excitations, a QED point-particle description again 
holds. For y/s ~ 2 m e , 2 m^, 2 m T the Z2 monopole is excited radially: A QED point- 
particle description fails. For y/s 3> mz the Hagedorn phase boundary of SU(2) e is 
locally overcome (the Z boson is interpreted as the decoupled dual gauge mode on 
the magnetic side of the magnetic-center phase boundary, the W bosons as the 
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decoupled gauge modes on the electric side of the electric-magnetic phase bound- 
ary): The multiplicity of final states should be in stark contradiction to the SM 
prediction (we expect charge nonconservation in such processes). Where does the 
parity violation come from? This is an intermediate consequence of the existence of a 
Planck-scale axion. What is the nature of the charge-zero state (single center- vortex 
loop)? The mass of this state for SU(2) e is roughly given by m v ~ compare 
with Eq. (|168p . Moreover, the mass of the Z boson is given as mz ~ gdec Tn e . From 
the experimentally known values m e ~ 5 x 10 5 eV and mz ~ 9 x 10 10 eV we thus 
have gdec ~ I x 10 5 and therefore m v ~ ^ eV. This is close to the upper bound 
for the mass of t he electron neutrino obtained from a tritium [3 decay experiment: 
m v < 2.3 eV ^^ l Thus a single center-vortex loop of SU(2) e viably is a candidate for 
the electron neutrino. Notice that this soliton has no antiparticle: Neutrinos need to 
be of the M ajor a na type in accord with the successful search for neutrinoless double 
beta decay A similar situation holds for SU(2) M and SU(2) T . We expect the 
masses of their intermediary vector bosons mz' ,m w >,± and mz", m W n,± to scale 
with their Yang-Mills scales m M and m T and large values of the gauge couplings 
at the respective phase boundaries. Thus there are very weak and very, very weak 
interactions in addition to the weak interactions which, however, will be very hard 
to detect experimentally. (To detect, say, mz 1 directly would need a center-of-mass 
energy in e + e~ annihilation which should at least be two- hundred times mz.) A 
remark concerning point (v) is in order: Since SU(2) e x SU(2) A1 x SU(2) r are in their 
center phases at present their contribution to the ground-state energy density and 
pressure of the Universe is nil. Above, we have computed the contribution arising 
from SU(2) cmb when assuming the Planck-scale axion to be absent. 

Let us make a short remark on the solar wind. This particle flux is mainly com- 
posed of protons (about 3 x 10 43 protons depart annually from the solar surface 
HQ2I). If the conservation of electric charge, which is a built-in feature of the SM, 
would hold then the sun would continuously acquire negative charge: A disastrous 
implication for earth's orbit would arise. The problem is resolved by the observation 
that in the solar core temperatures are greater than m e . Electronic charge, however, 
is absent in the magnetic or electric phase of SU(2) e . According to the phase dia- 
gram of SU(2) e the solar core contains a superconducting mantel (magnetic phase, 
Bose condensate of electric monopole-antimonopole pairs) whose negative pressure 
together with gravity balances the positive thermal pressure of the innermost core 
(electric phase) where fusion takes place. Within the core region there is a clear 
dominance of positive charge which the sun deposes off by m eans of the solar wind. 
A superconducting core of the sun is also demanded in for other reasons. 

We conclude this section by stressing that no fundamental Higgs field is needed 
to break the weak symmetry S\J(2)w or SU(2) e . In contrast to the SM, where this 
symmetry breaking is complete by a nonvanishing expectation of a fundamental 
Higgs field, the dynamical breakdown of SU(2) e proceeds in a two-stage, Higgs- 
particle free way: SU(2) e — >U(1) (electric phase; adjoint nonfluctuating Higgs field) 
and U(1)d — > 1 (magnetic phase; complex nonfluctuating Higgs field). Moreover, 
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the electron and its neutrino are stable solitons in the center phase of SU(2) e . 

7.3. Quantum chromodynamics : Quark confinement and fractional 
quantum Hall effect 

In this section we pursue an admittedly speculative and not very matured idea 
about the nature of quarks and their interactions. 

Quantum chromodynamics (QCD) is an integral part of the SM. QCD is the 
gauge theory of strong interactions: Pointlike current quarks, which are spin-1/2 
fermions of to-be-measured masses, are fundamentally charged under the gauge 
group SU(3)c and interact by the exchanges of massless gluons, the gauge bosons 
of SU(3)c- The latter interact with one another according to a pure Yang-Mills 
action. The electric charges of quarks are 2/3 or —1/3. Let us only discuss the three 
quark flavors of lowest mass m u = (3 • • • 5) MeV (charge 2/3), m<i — (5 • • • 7) MeV 
(charge -1/3), and m s = (100- • -140) ~MeV (charge-1/3) (all MS scheme, results 
depend on the renormalization point). 

Since leptons are likely to be the stable solitons in the center phases of pure SU(2) 
Yang-Mills theories it is tempting to speculate that quarks are related to the charge- 
one solitons (center- vortex loops with one self-intersection, spin-1/2 fermions) in the 
center phases of pure SU(3) Yang-Mills theories. If we assign an SU(3) U and SU(3)d 
theory of Yang- Mills scale ~ m u and rrid to the quark flavors u and d, respectively, 
then we need to understand why these quarks are confined and why the electric 
charge appears to be 2/3 or —1/3. A plausible way of generating quark confinement 
would be to add an additional SU(3) Yang-Mills theory of scale, say A = 140 MeV 
which, however, is a magnetic dual to the other SU(3) theories. A center-vortex 
condensate of this theory constrains the (color)electric flux between auord quark 
and a u or d antiquark into a tube and thus confines. The additional SU(3) theory 
also has charge-one solitons in its center phase which are confined by the center- 
vortex condensates of the theories SU(3) U and SU(3)<j. It is tempting to interpret 
these solitons as strange quarks s and thus to invoke the label SU(3) S . 

What about the electric quark charges? Due to confinement the trajectories of 
each quark flavor are forced onto a more or less two-dimensional spherical surface 
if the ground state of a given hadron is considered. A flux dual to the flux in the 
confining tube is readily available in the center-vortex condensate being responsible 
for confinement. This is the set-up for the occurrence of the fractional quantum Hall 
effect: Quarks that would be integer charged spin-1/2 fermions in the absence of 
the dual fluxes form bound states with these fluxes. Bound states with three dual 
flux quanta are bosons, and thus they condense. Excitations above this condensate 
have fractional electric charge. For a thorough discussion of this phenomenon, see 

EDS 

Again, all of what was said in this section is preliminary. We would like to stress 
though that the above scenario has the potential to explain why the equation of 
state in hydrodynamical simulations of the elliptic flow measured in ultra-relativistic 
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heavy-ion collisions at RHIC seems to be so close to the free- gas li mit despite the fact 
that strong correlations thermalize the system very rapidly lES] (At T c ~ 170 MeV 
the two theories SU(3) U and SU(3)d are deep inside their electric phases while 
SU(3) S is just above the electric-magnetic transition.) 

8. Conclusions 

We have developed a nonperturbative approach to SU(2) and SU(3) Yang-Mills 
thermodynamics. The formation of a macroscopic, adjoint, and nonfluctuating Higgs 
field in the deconfining (electric) phase of each theory, which involves the (admissible 
part of the) moduli space of a caloron-anticaloron system, the Bose condensation of 
thermalized magnetic monopole-antimonopole systems into a macroscopic, nonfluc- 
tuating complex field in a preconfining (magnetic) phase, and the nonthermal con- 
densation of systems composed of a center- vortex and an anti-center- vortex in a con- 
fining (center) phase were shown. A change of the statistics of the excitations from 
bosonic to fermionic was observed across the last phase transition. The degeneracy 
of the ground state with respect to a (global) electric Z 2 (SU(2)) and Z 3 (SU(3)) 
symmetry was observed in the electric phase, and the uniqueness of the ground state 
with respect to these symmetries was derived in the magnetic phase. Moreover, the 
nature of the phase transitions, electric-magnetic and magnetic-center, was clari- 
fied. The evolution of thermodynamical quantities with temperature was computed 
within the electric phase and the magnetic phase, and the density of fermionic states 
was estimated in the center phase. 

It did not escape the author's attention that the results obtained may resolve a 
number of long-standing problems in particle physics and cosmology. 
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